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SECTION  1 
INTRODUCTION 


This  manual  is  written  to  give  the  potential  user  of  the  DRAG  program  the 
specific  information  he  or  she  will  need  to: 

a)  Set  up  the  input  file 

b)  Run  the  program 

c)  Interpret  the  results. 

Prior  to  discussing  these  things  however,  let  us  present  some  background 
information  on  the  program. 

DRAG,  is  an  acronym  for  "Driver  Air  Cushion”  and,  as  the  program  name 

indicates,  the  program  was  written  to  describe  the  interaction,  during  a crash, 
between  the  driver  of  a vehicle  and  an  airbag.  Other  programs  have  been  written 
to  describe  such  an  interaction,  but  none  were  specifically  suitable  to  the  needs 
of  this  project  which  are: 

1)  To  have  the  capability  of  simulating  an  airbag  shape  typical  of  the 
ellipsoidal  shape  which  almost  all  driver  bags  have. 

2)  To  be  sufficiently  simple  and  inexpensive  that  it  can  run  on  a small 
computer  with  no  "library  routines"  necessary  for  execution. 
Therefore,  the  program  must  be  self  contained  and  inexpensive  to 
operate  so  it  can  be  used  as  a design  tool. 

3)  To  be,  as  a design  tool,  oriented  to  the  user  requirements  of  a typical 
restraint  system  engineer,  with  both  the  f emulation  and  the  input  - 
output  in  units  commonly  used  and  measured. 

4)  To  be,  as  a design  tool,  oriented  toward  the  test  hardware  actually 
encountered  in  most  situations.  For  example,  past  computer  programs 
might  model  the  driver  very  well  but  neglect  the  bag  shape  actually 
used  and/or  the  column  binding  and  frictional  forces  which  are  almost 
always  present  and  influence  the  results  greatly. 


1 


5)  To  have  the  same  level  of  detail,  complexity  and  accuracy  for  all  the 
components  of  the  restraint  system.  For  example,  there  are  several 
different  kinds  of  steering  colums,  all  of  which  behave  differently. 
It  does  little  good  to  have  an.  elaborate  airbag  algorithm  and  then  to 
back  it  up  with  a steering  column  that  is  so  simple  it  is  only 

described  by  a single  force- stroke  characteristic.  Such  a program 
must  at  least  have  the  capability  of  describing  the  different  column 
frictional  and  binding  effects  which  are  so  important  in  determining 
the  driver's  injury"  levels. 

In  writing  the  DRAG  computer  program  we  have  sought  to  fulfill  these 

needs . 
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SECTION  2 

PROGRAM  DESCRIPTION 


DRAT  is  a two-dimensional,  limped  mass  computer  model  of  the  driver 
interacting  with  an  ellipsoidal  airbag  mounted  to  a stroking  steering  column. 
The  model  includes  three  masses,  representing  the  head,  torso  and  lower  body. 
The  airbag  is  simulated  by  a ellipsoid  into  which  a programmed  amount  of  gas 
flows.  By  adjusting  the  airbag  vent  size,  a selected  amount  of  gas  can  be  vented 
during  the  crash  in  order  to  attenuate  the  peak  chest  Gs  and  the  rebound  effects. 

The  model  also  includes  the  steering  column  force-stroke  characteristics 
(including  all  binding  and  frictional  effects),  the  neck  rotational  resistance 
and  the  seat  friction. 


2.1  PROGRAM  FORMULATION 

The  mathematical  formulation  of  the  equations  of  motion  follows  the  classical 
Lagrangian  derivation  (Appendix  A of  this  manual)  with  body  pivot  points  at  A and 
B of  Figure  1.  The  lower  body  mass  (hips  and  legs)  is  constrained  to  move 
horizontally. 

DRAG  uses  a fixed  time  step  integration  routine  to  solve  the  differential 
equations  of  motion  numerically.  The  integration  routine  chosen  was  the  Adams- 
Moulton  predictor  corrector  method,  with  the  fourth  order  Runge-Kutta  method 
employed  to  determine  the  first  four  solution  points.  DRAG  has  been  written 
in  FORTRAN  IV  and  was  developed  to  operate  in  an  interactive  time  share  mode. 
The  program  is  self-contained,  in  that  no  external  routines  are  required  for 
execution.  It  is  also  modular  in  construction,  so  as  to  facilitate  the  addition 
of  other  subroutines  at  a future  time,  if  desired.  The  data  input  is  from  a 
previously  created  disk  file  and  the  input  parameters  appear  directly  on  the 
terminal,  immediately  preceding  the  complete  program  output. 
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FIGURE  1.  THE  MATHEMATICAL  MODEL 
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As  mentioned  above,  DRAG  has  been  programmed  in  modular  format  with  several 
subroutines.  The  reason  for  this  is  to  enable  other  features  to  be  easily  added 
at  a future  time.  For  example,  consider  the  tabular  data  input.  DRAG 
provides  two  basic  ways  to  operate  on  this  tabular  input.  For  simple  data  (where 
the  particular  value  of  the  dependent  variable  is  a function  only  of  the  value  of 
the  independent  variable),  a simple  table  look  up  and  interpolation  subroutine, 
"LOOKUP,"  is  provided.  Gas  flow  versus  time,  vehicle  acceleration  versus  time, 
neck  torque  versus  angle,  and  column  force  versus  crush  are  examples  of  this 
method  of  data  retrieval.  However,  in  those  cases  in  which  the  dependent 
variable  is  a function  not  only  of  the  value  of  the  independent  variable  but  also 
depends  on  whether  the  independent  variable  is  increasing  or  decreasing,  a 
different  subroutine,  "SPRING,"  that  allows  for  plastic  behavior,  is  used.  This 
subroutine  is  used  principally  for  those  cases  in  which  hysteresis  (or  plastic 
action  of  a deformable  member ) is  modeled.  In  this  case,  one  must  not  only 
specify  the  values  of  the  dependent  variable,  for  different  values  of  the 
independent  variable  but  must  also  specify  the  "unload  slope"  for  those 
conditions  in  which  the  member  is  undergoing  unloading  during  a lessening  of  the 
degree  of  deformation.  Knee  restraint  force  (as  a function  of  crush)  and  seat 
friction  (as  a function  of  stroke)  are  handled  by  this  subroutine. 


2.2  AIRBAG 

Most  of  the  driver  airbag  interaction  models  being  used  today  rely  on  relatively 
simple  spherical  or  cylindrical  airbag  shapes  (In  which  the  bags  exhibit  a 
constant  radius  of  curvature  regardless  of  impact  angle).  Unfortunately,  these 
simple  bag  shapes  do  not  adequately  describe  the  shape  of  most  of  the  driver 
airbags  presently  being  used. 

Of  all  the  geometric  shapes  that  could  be  postulated  as  candidates  for  the  driver 
airbag,  the  ellipsoid  is  most  nearly  the  shape  of  the  driver  airbag.  For  this 
reason,  we  chose  the  ellipsoid  as  the  shape  upon  which  to  base  the  bag  shape 
algorithm  described  in  Appendix  B of  this  manual.  Unfortunately,  the  ellipsoid 
is  not  as  mathematically  easy  to  describe  as  the  sphere  or  cylinder,  since  the 
bag  radius  of  curvature  and  the  intercepted  volume  of  the  airbag  are 
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Figure  2 shows  a simple  schematic  of  the  airbag,  along  with  the  variables 
necessary  to  describe  the  airbag  shape.  The  airbag  is  assumed  to  be  symmetric 
about  a line  coincident  with  the  steering  column's  longitudinal  axis. 


2.3  STfflIING  COLUM 

As  previously  mentioned,  the  proper  modeling  of  the  steering  column  is  very 
important  to  the  overall  accuracy  of  the  program  in  reproducing  ’’real  world” 
behavior  for  the  driver  restraint  system.  Most  driver  restraint  models  specify 
only  a force  versus  stroke  characteristic  for  the  colom.  In  many  cases  these 
models  neglect  column  mss,  column  frictional  and  binding  effects,  and  the 
specific  points  at  which  the  column  is  supported.  In  setting  up  the  column  for 
this  program,  we  have  chosen  a generic  type  of  colimn  that  is  widely  used.  This 
column  is  the  General  Motors  type  colum  shown  in  Figure  3.  Other  types  could 
have  been  specified,  but,  because  of  the  modular  program  construction,  other 
types  may  be  easily  added  as  required.  As  it  turns  out,  the  GM  type  column  is 
fairly  typical  of  a wide  variety  of  columns,  so  the  program  can  be  used  "as  is” 
for  a wide  variety  of  situations. 

The  program  is  set  up  to  calculate  the  airbag  forces  and  pressures  before  coming 
to  the  colum  force  calculation.  Therefore,  the  airbag  forces  are  resolved  into 
a column  axial  force,  a column  normal  force  and  a colum  moment  (as  applied  loads 
to  the  colum)  before  entering  the  column  dynamics  routine.  All  the  user  need 
specify  is  the  colum  mass,  the  column  angle,  the  basic  force-  stroke  properties 
of  the  column,  the  coefficient  of  friction  at  the  colum9 s support  points,  and 
the  pertinent  colum  dimensions  (as  shown  in  Figure  3)  for  the  computer  to 
calculate  the  complete  column  dynamics. 

Appendix  C of  this  manual  contains  the  details  of  the  steering  column  algorithm. 
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FIGURE  3.  04  ACRS  COLUMN 
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2.4  KNEE  RESTRAINT  AND  SEAT  FRICTION 


The  program  is  set  up  to  accept  tabular  input  for  the  force  versus  crush  - 
properties  of  the  knee  restraint  and  the  force  versus  displacement  properties  of 
seat  friction.  It  is  these  values  which  determine  what  the  lower  body  restraint 

for  the  lower  torso,  hips  and  legs  will  be. 

The  user  specifies,  in  tabular  format,  what  these  properties  are  to  be.  In 
addition,  the  user  specifies  the  "unload  slope,”  so  the  program  can  compute  the 
unloading  force  path  to  be  taken,  during  rebound,  away  from  the  knee  bolster  or 
movement  rearward  across  the  seat.  Specifics  of  how  the  input  is  handled  are 

discussed  in  Section  3.1. 


2.5  THE  DRIVER 

The  driver  is  modeled  by  three  lumped  masses  (the  head,  the  torso  and  the  lower 
body)  which  pivot  at  points  A and  B in  Figure  4.  Figure  4 also  describes  the 
driver  geometry  and  the  location  of  the  airbag  and  the  driver  with  respect  to  the 
compartment * Specific  details  needed  to  provide  driver  related  input  to  the 

program  are  described  in  Section  3.1. 


A eoraient  is  necessary  on  the  resisting  torque  generated  by  neck  muscular 
resistance  and  the  anatomical  interferences  caused  by  relative  displacement 
between  the  head  and  torso.  These  input  values  are  only  applied  if  certain 
conditions  are  met.  Thus,  in  cases  in  which  the  head  is  returning  to  be  more 
nearly  in  line  with  the  chest  (0^-0^  becoming  smaller)  the  torque  is  not  applied. 
Only  when  the  head  is  going  more  out  of  line  with  the  chest  (0^-Gtj,  becoming 
larger)  is  the  neck  torque  applied;  that  is, 


For  0H  - ©T  > 0 
For  0H  - ©T  > 0 
For  0||  - ©j  < 0 
For  0^  - ©«j.  < 0 

where: 


and  % ■ ©y  > 0>  T < 0 
and  ©u  • < 0,  T - 0 

and  ~ T s 0 

and  0H  - ©T  < 0,  T > 0 
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Neck 

Rotation 

Resistance 


T = Neck  resisting  torque 

©n  is  as  defined  in  Figure  4 

Oj  is  as  defined  in  Figure  4 

©H  - e.p  is  as  defined  in  Figure  1 

©n  - ©T  is  the  time  rate  of  change  of  ©H  - ©^ 

In  addition  to  the  neck  resisting  torque  (based  on  the  angular  displacement  of 
the  head  relative  to  the  torso),  we  also  have  included  a damping  coefficient  for 
the  head  (based  on  the  angular  velocity  of  the  head  relative  to  the  torso).  This 
value  is  known  as  "DCN*’  in  the  input  data.  Thus  the  overall  resistance  of  the 
head  to  rotation  with  respect  to  the  torso  is  composed  of  two  terms  — one  term 
based  on  relative  angular  displacement  and  a second  term  based  on  relative 
angular  velocity. 

Appendix  D of  this  manual  is  an  overall  listing  of  the  DRAG  program. 
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SECTION  3 

AN  EXAMPLE  OF  THE  PROGRAM'S  USE 


In  order  to  illustrate  the  use  of  the  program,  we  will  describe  a simulation  of  a 
sled  test  (Minicars  Test  1411)  of  the  Large  Research  Safety  Vehicle  (LRSV). 


3.1  CREATING  THE  DATA  INPUT  FILE 

The  first  thing  one  must  do  in  preparation  for  making  a computer  run  is  to  set  up 
a data  input  file  of  information  which  the  program  needs  to  run.  Table  1 is  a 
list  of  the  parameters  that  make  up  this  file. 

In  Table  1,  the  first  column  shows  the  location  of  a particular  piece  of  data  in 
the  input  file.  For  example,  1-1  refers  to  the  first  data  pointof  Line  1,  1-2 
refers  to  the  second  data  point  of  Line  2,  etc.  The  second  column  lists  the 
alphanumeric  name  of  the  variable  at  that  particular  place  in  the  data  file  and 
the  units  that  the  variable  must  have  in  the  input  file.  The  third  colum 
contains  a short  description  of  what  the  variable  is,  and  the  forth  column  gives 
the  actual  value  of  the  variable  for  sled  test  1411,  as  it  was  input  into  the 
data  file.  Thus,  Table  1 gives  all  the  required  input  to  simulate  test  1411. 
The  file  resulting  from  this  input  is  shown  in  Figure  5.  We  have  chosen  to  call 
the  file  LRSV. 
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TABLE  1.  INPUT  FILE 


Variable 

Location 


Line 

No. 

Location 
In  Line 

Variable  Name 
and  Units 

Description 

Value 

1 

1 

y(4)  mph 

Vehicle  impact  velocity 

39.3 

1 

2 

y(6)  degrees 

Head  angle  (0  in  Figure  1) 

-5 

I 

3 

y(7)  degrees 

Torso  angle  (0T  in  Figure  1) 

-18 

2 

1 

pounds 

Weight  o£  lower  body 

77 

2 

2 

Z^  pounds 

Weight  q£  torso 

67 

2 

3 

Z^  pounds 

Weight  of  head 

11 

2 

4 

Rj,  inches 

Distance  from  H-point  to 
torso  center  of  gravity 
(Figure  4) 

13.8 

2 

5 

inches 

Distance  from  H-point  to ' 
neck  pivot-  (Figure  4) 

19.5 

2 

6 

inches 

Distance  from  neck  pivot 

to  head  center  of  gravity 

when  0 = 0 (Figure  4) 
n I 

4 

2 

7 

rtoph  inches 

Distance  from  H-point  to 
top  of  head 

28 

3 

1 

NPN 

No.  of  points  on  Lines  14 
and  IS 

9 

3 

2 

NKR 

No.  of  points  on  Lines  18 
and  19 

8 

3 

3 

NV 

No.  of  points  on  Lines  16 
and  17 

6 

3 

4 

NSF 

No.  of  points  on  Lines  12 
and  13 

6 

3 

5 

NPG 

No.  of  points  on  Lines  4 
and  5 

8 

3 

6 

NPC 

No.  of  points  on  Lines  6 
aid  7 

6 

3 

7 

SUN  lb/ in 

Seat  friction  unload  scope 

1000 
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TABLE  1.  (CONT'D) 


Variable 

Location 


Line 

Location 

Variable  Name 

No. 

In  Line 

and  Units 

Description 

Value 

3 

8 

SKR  lb/ in 

Knee  restraint  unload  scope 

2900 

4 

1 

(234(1,1}  msec 

Gas  flow  time,  first  point 

Figure  5, 

4 

etc. 

GEN(l,k)  msec 

Gas  flow  time,  subsequent 

Line  4 

points 

5 

1 

GEN(2,1)  lb/ sec 

Gas  flow  rate,  first  point 

Figure  5, 

5 

etc. 

GEN(2,k)  lb/sec 

Gas  flow  rate,  subsequent 
points 

Line  5 

6 

1 

C0L(1,1)  in 

Column  stroke,  first  point 

, Figure  5 , 

6 

etc. 

COL(l,k)  in 

Column  stroke,  subsequent 
points 

Line  6 

7 

1 

C0L(2,1)  lb 

Column  force,  first  point 

Figure  5, 

7 

etc. 

C0L(2,k)  lb 

Column  force,  subsequent 
points 

Line  7 

8 

1 

ATMOP  PSIA 

Local  atmospheric  pressure 

14.7 

8 

2 

PGZ  PSI G 

Initial  airbag  pressure 

-14.7 

8 

3 

GTZ  Deg  R 

Temperature  of  gas  entering 
airbag 

1160 

8 

4 

u ib  \ 

Gas  constant 

660 

8 

5 

m 

PN1 

Polytropic  gas  exponent, 
flow 

1.4 

8 

6 

PN2 

Polytropic  gas  exponent, 
compression 

1.4 

8 

7 

PN3 

Polytropic  gas  exponent, 
expansion 

1.4 

9 

1 

va 

Vent  discharge  coefficient 
- subsonic  flow 

0.7 

9 

2 

VC2 

Vent  discharge  coefficient 
- sonic  flow 

0.7 

9 

3 

AV  sq  in 

Vent  area 

3.5 

9 

4 

SA  inches 

Major  axis  length  of 

11.5 

airbag  (Figure  2) 
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TABLE  1.  (CONT’D) 


Variable 

Location 

Line 

Location 

Variable  Name 

No. 

In  Line 

and  Units 

Description 

Value 

9 

5 

SC  inches 

Minor  axis  length  of 
airbag  (Figure  2) 

7.5 

9 

6 

X1Z  inches 

Horizontal  rim  reference 

dim.  (Figure  4) 

29.5 

9 

7 

YlZ  inches 

Vertical  rim  reference 

dim.  (Figure  4) 

23.0 

9 

8 

ft-lb-sec 

Neck  damping  coefficient 

2.5 

DLJM  -i 

rad 

10 

1 

SLIM  inches 

Column  stroke  limit 

8.0 

10 

2 

THETAC  degrees 

Column  angle  (Figure  3) 

17 

10 

3 

MU 

Column  frictional  co- 
efficient (Figure  3) 

0.16 

10 

4 

LSCZ  inches 

Column  reference  dirft.  to 
shear  capsule  (Figure  3) 

16 

10 

5 

LFWZ  inches 

Column  reference  dim.  to 
firewall  (Figure  3) 

14.75 

10 

6 

LBAZ  inches 

Column  reference  dim.  to 
aft  bushing  (Figure  3) 

16 

10 

7 

LBFZ  inches 

Column  reference  dim  . to 
forward  bushing  (Figure  3) 

23.25 

10 

8 

WC  pounds 

Column  stroking  weight 

18 

11 

1 

WH  inches 

Head  width 

9 

11 

2 

YO  inches 

Y'  coordinate  point  "0" 
(Figure  2)  always  = 0 

0 

11 

3 

RIMRAD  inches 

Radius  of  steering  wheel 
(Figure  3) 

7.75 

11 

4 

X2Z  inches 

Horizontal  H-point  reference 
dim.  (Figure  4) 

36.5 

11 

5 

Y2Z  inches 

Vertical  H-point  reference 
dim.  (Figure  4) 

9 

11 

6 

WB  inches 

Width  of  torso 

17 

13 


TABLE  1.  (CQNTB) 


Variable 

Location 


Line 

No* 

Location 
In  Line 

Variable  Name 
and  Units 

Description 

12 

1 

SFM(1,1)  inches 

] 

Seat  friction  displacement, 
first  point 

12 

etc. 

SFN(l,k)  inches 

Seat  friction  displacement, 
subsequent  points 

13 

1 

SFN(2,1)  lb 

Seat  friction  force,  first 
point 

13 

etc. 

SFN(2,k)  lb 

Seat  friction  force, 
subsequent  points 

14 

1 

FNECKU.l)  deg. 

Head  relative  angle  (0^  -0^, 
Figure  1) , first  point 

► 

14 

etc. 

FNECK(l,k)  deg. 

Head  relative  angle, 
subsequent  points 

15 

1 

FNECK(2,1)  ft-lb 

Neck  rotational  resistance 
torque,  first  point 

15 

etc. 

FNECK(2,k)  ft-lb 

Neck  rotational  resistance 
torque,  subsequent  points 

16 

1 

VEHGS(ljl)  msec 

Crash  pulse  time,  first  point 

16 

etc. 

VEHGS(1 ,k)  msec 

Crash  pulse  time,  subsequent 
points 

17 

1 

VEHGS(2,1)  Gs 

Crash  pulse  Gs,  first  point  1 

17 

etc. 

VEHGS(2,k)  Gs 

Crash  pulse  Gs,  subsequent 
points  J 

18 

1 

KRN(1,1)  in 

Knee  displacement,  first 
point 

18 

etc. 

KRN(l,k)  in 

Knee  displacement, 
subsequent  points 

19 

1 

KRN(2,1)  lb 

Knee  (femur)  force  total, 
2 knees,  first  point 

19 

etc. 

K3RN(2,k)  lb 

Knee  (femur)  force  total, 

2 knees,  subsequent  points 

Value 


Figure  5, 
Line  12 


Figure  5, 
Line  13 


Figure  5, 
Line  14 


Figure  5, 
Line  15 


Figure  5, 
Line  16 


Figure  5, 
Line  17 


Figure  5, 
Line  18 


Figure  5, 
Line  19 
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FIGURE  5.  TOE  INHJT  FILE  LRSV 


3.2  RUNNING  TOE  PROGRAM 

Once  the  input  file  has  been  created  and  saved,  we  are  ready  to  run  DRAG  • At 
this  point  the  user  accesses  DRAG  and  tells  it  to  run.  The  computer  will 
respond  by  asking  the  user  to  name  the  input  file;  in  this  case  we  respond 
"LRSV,”  as  shown  in  Appendix  E of  this  manual.  Once  this  answer  is  given,  the 
program  will  begin  to  print  out  the  input  data  - first  in  the  units  input  into 
the  program  and  then  in  the  units  used  by  the  program  in  the  actual  computation. 

Let  us  now  discuss  the  output. 

Altogether  there  are  seven  blocks  of  output  — each  block  consisting  of  the 
amount  of  data  that  can  be  conveniently  grouped  together  in  terms  of  subject 
Block  1-of  output  in , App . E has  Elapsed  Time,  Chest  A-P  Acceleration, 
Chest  S-I  Acceleration,  Head  A-P  Acceleration  and  Head  S-I  Acceleration. 


IS 


Block  2.-of  output  in  Appendix  E is  basically  in  sumary  of  the  main 


parameters  of  interest  in  the  program  and  consists  of  the  following  items 
(page  78  )t  Elapsed  Time?  Vehicle  Gs  (Acceleration) , Vehicle  Velocity, 

Vehicle  Displacement  (Crush) , Body  Gs  (chest  A-P  Gs) , Column  Displacement 
(Stroke) , and  Airbag  Pressure. 

Block  3 - (page  79  ) has  Elapsed  Time,  H-Point  Displacement  (with  respect  to 
ground),  H-Point  Velocity,  H-Point  Acceleration,  Femur  Force  (each  femur  has 
this  force).  Seat  Friction  Force,  and  H-Point  Relative  Displacement  (with 
respect  to  original  position  in  compartment) . 

Block  4 - (page  80  ) has  Elapsed  Time,  Torso  Displacement  (with  respect  to 
ground) , Torso  Angle  (with  respect  to  a vertical  line,  positive  when  toward 
dash) , TOrs © Angular  Velocity  (with  respect  to  H-point  pivot) , Torso  Angular 
Acceleration,  Torso  Relative  Displacement  and  Torso  Relative  Velocity. 

Block  5 - (page  81  ) is  the  exact  equivalent  of  Block  3,  but  for  the  head 
with  respect  to  the  neck  pivot  (instead  of  the  torso  with  respect  to  the 

H-point  pivot) . 

Block6  -(page  82  ) has  Elapsed  Time,  Column  Axial  Applied  Force,  Column 

Itoonal  Applied  Force  (positive  when  upward) , Column  Applied  Moment  (positive 
when  bending  the  column  upward) , Force  Resisting  Column  Stroke,  Column  Stroke, 

and  Column  Stroking  Velocity  (with  respect  to  the  compartment) . 

Block  7 - (page  8 3 ) has  the  Elapsed  Time,  Bag  Penetration  (measured  normal 
to  the  torso,  halfway  between  the  two  bag  intercept  points) , Bag  Volume,  Bag 
Pressure,  Bag  Wraparound  (fabric  tension)  Force  (on  the  chest)  and  Bag 

Pressure  Force  (on  the  chest) . 
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3.3  DRAG  OUTPUT  CORRELATION  WITH  TEST  RESULTS 


In  Section  3.1  we  described  the  use  of  data  for  an  LRSV  sled  test  as  input  to  - 
DRAG,  The  correlation  between  the  data  resulting  from  the  sled  test  and  the 
output  data  from  DRAG  is,  in  general,  a measure  of  the  predictive  capability 
of  the  simulation  model. 

Figures  6 through  10  illustrate  the  comparison  between  the  experimentally 
determined  test  results  and  the  DRAG  predictions  with  identical  input 
parameters . 

Legends  Solid  - test  data 


FIGURE  6.  DRIVER  HEAD  ACCELERATION 
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VERTICAL  (6s)  LONGITUDINAL  (Gs) 


FIGURE  7.  DRIVER  CHEST  ACCELERATION 
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RIGHT  FEMUR 
FORCE  (lbs) 


TIME  (msec) 

FIGURE  8.  AIRBAG  PRESSURE 


TIME  (msec) 
FIGURE  9.  FEMUR  LOADS 
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FIGURE  10.  COLUMN  COLLAPSE 


The  actual  HIC  for  the  sled  test  was  130  with  t^  = 42  msec  and 
=163  msec.  The  computer  simulation  calculated  an  HIC  of  158  with 
t^=50  msec  and  t2=130  msec. 
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APPENDIX  A 


DERIVATION  OF  THE  EQUATIONS  OF  MOTION 

The  derivation  of  the  equations  of  motion  will  be  formulated  utilizing  Lagran- 
gian  techniques  based  upon  the  geometrical  representation  in  Figure  A-l. 

Writing  an  expression  for  the  total  kinetic  energy  of  the  occupant,  we  have: 

(1)  T - ■>  0^  (%  + Y*)  + Mj.  (X^.  + Yp  + M^] 

Note  that  Yp  e 0,  as  no  movement  normal  to  the  X-direction  is  allowed  for  the 
hip -leg  mass. 

Mh  = Head  mass 

Mp  = Torso  mass 

Mp  = Hip -leg  mass 

Xp  = Horizontal  translation  of  the  hip -leg  mass  with  respect  to  inertial 
reference  point  — which  is  positive  when  it  is  in  direction  shown. 

• Xp  and  X^  are  similarly  defined 

Yjp  = Vertical  distance  from  H -point  to  the  center  of  gravity  of  the  head 
Yp  = Vertical  distance  from  H-point  to  the  center  of  gravity  of  the  torso 
Successive  dots  indicate  velocity  and  acceleration,  respectively. 

Writing  the  transformation  equations,  we  have: 


(2  J 

V 

*L  + rT  SIN9I 

(3) 

ii 

rT  coseT 

(4) 

V 

\ + rN  SIN6t  ♦ rH  SIN6h 

(5) 

V 

rN  C0S6T  + rH  C0S9H 

(6) 

V 

*L  + "I  C0S6T  ®T 

(7) 

V 

-rT  SLN6T  6t 
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MATHEMATICAL  MODEL 


Figure  A-l.  Geometrical  Representation  of  the 
Driver-Airbag  Interaction 
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(8) 


*H  = *L  rN  C0S9T  °T  + rH  C0S9H  ®H 

(9)  Yh  - -rN  SIN0t9t  - rH  SIN6H  0H 

where: 

r™  = Distance  from  hip  H-point  to  torso  center  of  gravity  (The  H-point 
is  assumed  to  be  coincident  with  the  hip -leg  center  of  gravity.) 

rN  = Distance  from  H-point  to  neck  pivot  point 

r„  = Distance  from  the  neck  pivot  point  to  the  center  of  gravity  of  the 
M head 

0H  and  are  as  defined  in  Figure  A-l. 

Substituting  Equations  6 through  9 into  Equation  1,  we  have: 

(10)  T - * {m,**  ♦ Hj  [x£  ♦ 2 XjTt  COS0t  0t  + 4 8*]  + ^ + 

2 XL  (rN  COS0T  9t  + rH  COS0H  + 2 rN  rH  (COS0T  COS0H  0T  0„  + 

SIN0t  SIN0h  9?  eH)  + ij  8*]J 

The  potential  energy  portion  of  the  Lagrangian  is: 

(11)  VT  - Mj.  g rT  COS0T 

C12)  VH  = ^ g (rH  COS9h  + rN  COS0T) 

Note:  The  applied  forces  and  moments  will  be  treated  separately  later  on. 

Writing  the  Lagrangian,  we  have: 

(13)  ' L = T-V  = T - (VT  + VH) , 

where  the  values  to  be  substituted  into  this  equation  are  given  by  Equations  10, 
11  and  12. 
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The  basic  equation  in  Lagrangian  mechanics  is: 

(14)  _d  / 3Lv  3L 

dt  \ 3q7"  / 3q^  qi 


where : 

vt.1L 

q^  s generalized  displacement  of  the  x mass 

th 

^ - generalized  velocity  of  the  i01  mass 

lr];l 

F . - generalized  force  acting  on  the  im  mass 

4X 

Taking  the  required  derivatives  from  Equation  13  for  substitution  into  Equa- 
tion 14,  we  obtain: 


(15)  3L 


a^  + + V \ + Nlr  rT  COS9t  9t  + ^ (rN  COS0t0t  + rH  C0S6h  6h> 


(16)  d_  / 3L_\  + Mj.  + jy  \ SINBt  9t  - (rN  SIN6T  9f  + 

dt  \ 3^  / 

rH  SIN6h  9^)  + Hj.  rT  COSeT  0T  + ^ (rN  COS0T  0T  + 


rH  coseH  eH) 


(17) 


(18) 


3L 

3Xl 

3L 

a§_ 


o 


Mp  (\  rT  coseT  + if  eT)  ♦ ^ 


\ rN  C0SeT  + rN  rH  ^C0SeT  C0SeH  9H 


+ SIN0T  SIN0H  eH)  + r*  e 


M 
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('19)  dt  ( ) = ^ rT  COS0T  ’ rT  SIN6T  9T  + rT  9T 


(21) 


(22) 


(23) 


+ ^ [ *L  rN  C0S6T  * K "N  SIN6I  9T  - rN  rH  (SIN6T  9T  C°S6H  9H 

+ coseT  sinqh  q2  - cos0T  cos0H  ©H  - cos0T  ©T  • sin©h  ©h 


- SIN0t  COS0h  02  - SIN0t  SIN0h  0h)  + r2  02 


] 


(20)  3L 


30  XL  rT  SIN9T  eT  ' rN  ^SIN9T  eT  XL  + rH  SIN0t  COS0h 

T 

9T  9H  “ rH  COS0T  SIN6H  9T  V + g rT  SIN9T  + ^ g rN  SIN9T 


3L 

30, 


= MH  [*LrHC0S9H 
+ rH  ®H] 

dt  ( 35^)  = “h  [\  rH  COS0H  " *L  rH  SIN0H 


+ rN  rH  (CC60t  COS6h  0t  + SIN0  SIN0  0„) 


®H  ‘ rN  rH  (SIN9T  C0S6H  0T 


+ cos0T  sin9h  eT  eH  - coseT  cos0H  0T  - cos©  sin©  e* 


SIN0t  COS0h  0t0h  - SIN0t  SIN0h  0t)  + r*  0^ 


- -^rH  rSIN0H9H\  + rN 
39h  L 

9T  V + ^ 8 rH  SIN6h] 


(COS0T  SIN0h  0t  Sh  - SIN0t  COS0h  X 


Substituting  Equations  16  and  17  into  Equation  14,  we  have: 

(24)  + Mj.  + Mjj)  - Mj.  rT  SIN0t  e|  - (rN  SIN0Te*  + rH  SIN©H  9*) 

+ Hf  rT  COS9T  0T  ♦ ^ (rN  COS0T  9'T  ♦ rH  COS6H  ^ 


which  is  the  equation  of  motion  for  mass  M^. 
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Substituting  Equations  19  and  20  into  Equation  14,  we  have: 


Hr  (XL  rT  COS0T  - *L  rT  SIN6T  0T  + rT2  V + “h  [ *L  TN  C0S8T 

• *L  rN  SINSI  0T  - rN  rH  (SIN6T  §T  C0S6H  ®H  + C0S6I  SIN9H 
• • 

- cos9t  coseH  eH  - coseT  eT  siNeH  eH  - sinst  coseH  9* 

- SIN0x  SIN0h  V + rN  0t]  + Hr  \ rt  SIN0t  0t  + “h  rN  ( 

SIN6t  9t  \ + rH  SIN6t  CQS©h  0t  ©h  - rR  COS0T  SIN0H  ©T  ©H) 

" Mj.  g rT  SIN0t  -fflgrN  SIN0T  = FQT  . 

Rewriting  the  above  yields: 

(25)  Mj.  (XL  rT  COS9t  + r*  9T)  + 1^  ^ rN  COS8T  - rN  rH  (COS8T  SIN8H  0* 

• • •• 

- coseT  coseH  a - sin8t  cosaH  e*  - siN8  sin8  e ) 


rN  9T 


Mj.  g rT  SIN9t  - Mjj  g rN  SIN8T  = F8. 


whidi  is  the  equation  of  motion  of  the  torso  mass. 

Substituting  Equations  22  and  23  into  Equation  14,  we  have: 

^ [Vh  COS0H  - *L  rH  SIN0H  9H  - rN  rH  (3IN9I  COS0H  0T 

• • 

+ C0S8t  SIN8h  8t  9h  - C0S8T  C0S8H  0T  - COSGj,  SIN9H  6* 

• e ••  ~ 

- SIN0J,  COS9h  8t  8h  - SIN9t  SIN8h  9?)  + r*  9H  + 1^  rH 

l « 

SIN0H  0h  ^ + rN  (COS0T  SIN0h  ©t  ©h  - SIN0T  COS0H  0T  0H) 


Mn  g rH  SIN0h  = F 


0H 
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Rewriting  the  preceding  yields: 


(26)  [xLrH  COS0H  - rN  rH  (SIN6T  C0S6H  6*  - COSeT  COS0H  0T 

- COS0T  SIN0h  0*  - SINeT  SIN0h  0t)  + rH  eH  1 

- ^ g rH  SIN6h  = F0h  , 

which  is  the  equation  of  motion  for  the  head  mass. 


Writing  Equation  24  in  terms  of  XL,  we  have: 

(27)  = 1 fFXL  + rT  + ^ V SINeT  + % rH  ® 

L Mj,  + ^ l 

- (Mr  rT  + Vn3  9t  C0S9t  - % rH  9h  coseH  } • 


Writing  Equation  25  in  terms  of  0^.,  we  have: 

0T  " Mj,  r*  + Mjr  r*  (F0T  " ^ rT  + “h  V *L  COS0T 

' ”h  rN  TH  [ SH  (COS0H  COS0T  + SIN0H  SIN9T) 

+ (-SIN6H  COS0t  + COS0H  SIN0t)1 
+ Hj.  g rT  SIN0t  t^grN  SIN0t  l . 
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cm  DC 


Writing  Equation  26  in  terms  o£  0^?  we  have: 


(29) 


9h  = 


0H 


-XL 


COS0 


Vii 


- ' — r. 


•H 


‘H 


(COS0U  COS0,,,  + SIN0  SIN0  ) 0 


H 


+ (sin0h  cos0t  - coseH  sin0t)  0"  1 + smeH 

j rH 


where: 


H 


F6H  ■ FH  rH  + T NECK 

F0T  = FH  C0S(eH  ‘ 6T)  rN  + FT  rF  • T “** 
FXL  = FH  C0S6H  + FT  C0SeT  + FL  * 
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APPENDIX  B 


DERIVATION  OF  THE  AIRBAG  ALGORITHM 


The  angles  and  distances  described  in  this  appendix  will  be  those  depicted  in 
Figure  B-l. 


COORDINATE  TRANSFORMATION  EQUATIONS 


(1) 

y = y2  + y" 

(2) 

xsx^  + x" 

(3) 

y = yx  + y' 

COS0  + X’ 
c 

SIN0 

t 

(4) 

X = x^  + x' 

COS0  - y* 
c 

SIN0 

t 

y 


To  obtain  transformation  equations  for  x"  aid  yIT  into  the  xr  ,y*  system, 
substitute  Equations  1 and  2 into  Equations  3 and  4 to  get: 


(5) 


xf 


x - x,  + x"  + y'  SIN0 
Z 1 c 

cose 

c 


(6) 


y' 


y - y + yM  - x’  SIN0 

w 1 C 

rase 

c 


Assume  that  the  torso  may  be  represented  by  a plane  that  intersects  the  airbag 
at  line  A-B  on  the  plane  of  symmetry  of  the  airbag  (as  shown  in  Figure  B-l) . 
Assume  further  that  the  airbag  is  an  ellipsoid  whose  plane  of  symmetry  in  the 
X-Y  plane  is  as  shown  in  the  figure.  Our  job  now  will  be  to  derive  an  equation 
for  the  bag  intercept  points  in  the  xT  -y r coordinate  system. 

In  the  x"-y"  system  the  equation  for  line  A-B  is: 

(7)  y"  = mx"  + b 


Equation  7 into  Equation  6 yields: 


y - y + mxM  + b - x'  SIN0 
2 1 c 

COS0 

c 


Substituting 

y'  ■ 
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COMPMTMEffr.  BM,  OilVER  COORDINATE  SYSTEM 


Figure  B-l.  Geometrical  Representation  o£  the  Airbag  Parameters 


Substituting  in  Equation  8 for  x"  from  Equations  2 and  4: 

- y in  (x'  COS0  - x_  + x_  - yf  SIN0  ) - x'  SIN0  + b 
. ‘Z  '1  c 2 1 c c 

(9)  y = cose ’ 


which  is  the  desired  equation  in  the  x'  ,yf  system. 

Let  y2  - y^  + b - m (x^  - x^)  = B (a  constant)  and  solve  Equation  9 for  y* 

y*  (COS0  + m SIN0  ) = B + x'  (m  COS0  - SING  ) 
c c c c 

B + x’  (m  COS0  - SIN0  ) 

(10)  y*  = COS0  + m SIN9~ ~ Equation  for  A-B  in  x\y'  system) 

c c 

The  equation  for  the  airbag  in  the  x'-y'  system  is: 


(11) 


I 2 yt2 
— + “ — = 1 


2 2 
c a 


Substituting  Equations  IQ  and  11  and  collecting  terms : 

xT 2 [a2  (COS0  + m SIN0  )2  + c2(-SXN0  + m COSO  )2] 
c c c c 

+ 2 Bc2x'  (m  COS0  - SIN0  ) + B2  c2  - a2  c2  (COSO  + m SIN0  )2  = 0 , 

c c c c 

which  is  a quadratic  equation  in  terms  of  x ' . 

2 2 2 2 
Let  A = a (COSO  + m SINO  ) + c -(SIN 0 t m COSO  ) 

C C 0 0 

D = 2 B C2  (m  COS0  - SIN0  ) 

c c 

E = B2  c2  - a2  c2  (COSO  + m SINO  )2 

c c 

Ax' ^ + Dx'  + E = 0 


(12)  x*  = 


-E>  + y D - 4AE 
2A 
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Values  for  x*  obtained  with  (12)  when  substituted  into  (10)  will  give  the  corres- 
ponding values  for  yf.  We  now  have  defined  the  line  of  intercept  (A-B)  of  the 
occupant's  body  with  the  mid-plane  of  the  airbag. 

With  this  line  now  established,  we  can  begin  to  calculate  the  restraint  forces 
that  will  be  applied  to  the  driver. 

Forces  will  now  be  calculated  due  to  pressure  effects  (Figure  B-2) . The  force 
on  the  head  aid  chest  are  composed  of  two  coup  orients  - a pressure  component  and 
a "wrap-around"  component  due  to  fabric  tension;  i.e.. 


(13)  F * F + F 
1 J QffiST  P FT 

c c 


(14)  F = F_  + F_ 

J HEAD  P„  FT 
ri  n 


Hie  pressure  forces  act  normal  to  the  head  and  chest: 


(15)  fp  -Pwb  <V  W 

c 


(rn  - W < “ 


P vl  AB 


(16)  Fp  "PwH 
H 


«ran  ' W < " 


pwh  [AB  ■ (rn  • wi 


«ioph  - W r " 


where  the  pressure  P must  be  calculated  due  to  bag  volume  and  thermodynamic 
effects. 

The  fabric  tension  component  will  be  calculated  later.  Let  us  now  calculate 
the  body  moments  caused  by  these  forces.  Using  the  H-point  and  neck  pivots  as 
our  reference  points: 


F0t  = FOiEST  * *FT 

(!8)  F - Fhead  * ^heaD  * 

H 


where  F^^  and  F are  given  by  Equations  13  aid  14.  We  will  now  evaluate 
Qihol  HEAD 

rft  rhead* 


32 


I 


I 


Figure  B-2.  Head-Bag  Interaction  Forces  Diagram 
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In  order  to  solve  for  R__,  we  must  derive  equations  for  the  H- point  location  in 
terms  of  the  xf-y’  coordinate  system  (see  Figure  B-3) . 


From  the  geometry  of  the  mid-plane  of  bag  impact,  the  H-point  coordinates 

are: 

(19)  x8h  - (y2  - yx)  SIN©c  + (x£  - xj  COS©c 

(20)  y*H  - C0S©c  [y2  - y±  - (x2  - x^  TAN©^]  . 

Ihe  equation  for  is: 


(21) 

*FT  ~ 

/(X,FT  • X'h)2  + ' 

where: 

(22). 

-V  9 _ 

(x’a  + x'neck-1 

i . . ^ 

FT 

2 

(23) 

V®  s 

7 FT 

y 5 + y * 

y A y NECK 
2 

The  equation  for  is: 

(24) 

rhead 

^TOPH  ‘ \ 
2 

(25) 

^HEAD 

2 

" + “mg  > % 

AB  + RMGiRN  ’ 


This  derivation  completes  the  solution  for  terms  needed  for  pressure  force  and 
body  moment  computation.  We  must  now  derive  equations  for  the  fabric  tension 
component  of  bag  force  due  to  bag  wraparound  in  the  lateral  plane.  (No  wrap- 
around in  the  vertical  plane  is  considered,  since  the  body  is  generally  as  long 
as  the  bag  is  high,  so  n©  wrap-around  will  occur.) 


34 


(y1-yo+(xo-x1)tan0r)s1ne 


T» 

s 

u» 
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Figure  B-3.  H-Point  Location  in  x'-y*  Coordinate  System 


Hie  figure  on  the  right  is  a 
view  looking  down  at  the 
deformed  bag. 


Let  us  now  consider  the  body  wraparound  forces  caused  by  fabric  tension.  This 
courponent  of  force  is  influenced  most  by  bag  pressure,  bag  penetration  and 

body  to  bag  width. 


Hie  force  N is  the  tensile  force  in  the  bag. 

<t> 

wb 

At  z 's  r—  , N is  obtained  by  a force  balance. 

2 <p 

N x perimeter  of  Section  A- A = P x Area  A-A , 
<P 

or 


(26)  N - 

<P 

Solve  for  a 8 


P x Area  A-A 
Perimeter  A-A 


w. 


and  c9  @ z 


(force  per  unit 
length  of  AB) 


We  may  write  for  the  y*,z9  plane, 

(27)  yj2  + z#2  * a2 

_ Wb 

For  z - Y 

I w 2 

(28)  y?  = a9  ~ ^ a2  - ~— 


36 


For  the  x,z  plane: 


Assume  the  area  of  Section  A-A  varies  linearly  with  bag  penetration  from 
its  initial  value,  ira’c*,  to  zero  when  fully  compressed  at  a bag  penetration 
of  G T.  Then: 


C3°)  Aa_a  = ira'C  (1  - ^ ) , 

where  BP  = bag  penetration  perpendicular  to  and  at  mid  point  of  torso 
GI  = length  across  bag  in  BP  direction. 

The  perimeter  is  given  by: 


(31) 


PER.  . ^ 2ir 
A-A  — 


(an  approximate  formula  with  accuracy 
sufficiently  close  to  the  exact  formula 
which  involves  elliptic  integrals) 


Substituting  (30)  and  (31)  into  (26) 


(32)  N = 


a’c'P 


J a1^  + c' 


Cl  --) 

U GI  J 


Substituting  (28)  and  (29)  into  (32) : 

T4 


P a 


2 Wb 


c2U 


(33)  N = 


_ w, 

.2  b 2 

2 ./  a — + c 

4 


B) 


r ( 


1 - 


BP 

GI 
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which  is  the  equation  for  the  tension  force  in  the  ellipsoidal  airbag  in  terms 
of  the  bag  pressure,  the  lengths  of  the  major  and  minor  axes,  the  bag  penetra- 
tion, and  the  driver  body  widths; in  units  of  force  per  unit  length  M. 


The  bag  perimeter  in  the  plane  at  the  mid- 
point of  AB  at  an  angle  u to  the  x-axis  can 
be  found;  see  figure  at  right.  This 
perimeter  will  remain  constant  and  the  wrap- 
around configuration  must  maintain  this 
perimeter.  To  find  the  perimeter  of 

ellipse  cut  by  B-B,  one  must  follow  the 
steps : 

1)  Find  the  distance  ST  (which  will 
be  the  length  of  one  axis) 

2)  Find  the  mid-point  of  GT 

3)  Find  the  other  axis  length. 


Calculate  BP  and  GI 
To  derive  the  equation  for  UT: 


(34)  Slope  of  GI  = • 

The  point  it  goes  through  is  Xp^,  ypj,. 
Then,  writing  the  equation  for  GT, 

(35)  (y*  - 7^.)  - ' XpP  ’ 


or,  rewriting  the  equation. 


(36)  y*  ~ + >Vr  * 


The  equation  for  the  mid-plane  of  the 
ellipse  is 


(37) 


Substituting  y!  into  the  above  and  collecting  terms , one  obtains 

^AB2  \ 2 ^AB  X’  ( , \ ,(yFT  “ "’PAB  ^Fl) 

2 ( yFT  " "VaB  *¥!  / 


(38) 


1 = 0 
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Let 


(39)  A1  = — y 


1 „ m'pAB2 

7 — T~ 

c a 


2 m 


B1 


PAS  , , . 

2 “ i/p T " m PAB  XFF 

2 


Cl  = ^yFT  " m>PAB  XFT^  _ x 


(40)  xG  - 


- V 


■B1  + > B1  - 4A1C1 

"ZA^ 


yG  ~ ^PAB  flXG  " XFT^  + yFT 


(41) 

Y S 

-B1  - M Br  - 4A1C1 

XI 

2^ 

(42) 

GI  = 

yI  " yG^  + ^XI  " XG^ 

(43) 

BP  = 

V (yFT  " yG^  + ^XFT  ” XG^ 

At  the 

midpoint  of  (ST, 

yI  ~ ^PAB  ^XI  " XFT^  + yFT 


*G  + XI 


*MGI 


y . V121 

yMGI  2 
In  the  plane  x*  yf , 

,2  .2 

S-V1  ; 

c a 

for  x*  = we  have: 

r2  = a2  ( 1 - ^ 


{'■4)  • 
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Now*  in  the  plane  y-z*  away  from  the  y?,  zf  axis,, 


2 

2 . 2 2 f - ^4GI 

z + y s a I 1 - — *— 


or»  for  y s yMc P 


2 

^MGI 


2 / yMGI  * 


which  is  the  length  of  the  other  axis* 


(45)  PerB  B - 2tt 


V 


— „ 2 2 
CGI/2)  + z 


where  GI  and  z are  given  by  Equations  42  and  44,  respectively. 

For  constant  perimeter  for  section  B-B, 


„ . RIM 

SBR  s 2 I 1^)  + DBR" 


PerBm  - w,  - SBR 
(46)  - V 


chord  length  at  y^;  RIM  * 2 ^ - yj  , where  R -a^  — .) 
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(47)  2R,SIN(180  - <f>/2)  = 2R  SINcf)/2  = 

<P  <P 


GI  - BP  - DBR 
COSa 


(48)  a = TAN 


»b 


- RIM 


2(GX  - BP  - DBR) 

Solving  (46)  § (47)  simultaneously  gives: 

(49) 

T ‘ -1 / 

ms  tam  b 

k2(GI  - BP  - DBR)^ 

Equation  (49)  must  be  solved  numerically  for  <f>. 

The  fabric  tension  force  can  now  be  calculated. 

(50)  F_,  = 2N  AB  COS (180  - (4>/2  +a)) 

FT  < p 

= -2NAB  COS  (<p/  2 +a) 

<P 

= -2NAB  C0S6, 

<P 


where  is  given  by  Equation  33.  AB  is  given  by  AB  - 


+ CyA  ' V2 


<f>  is  given  by  Equation  49,  and  a is  given  by  Equation  48. 

This  completes  the  restraint  force  computation.  It  remains  to  compute  the 
volume  change  as  a function  of  bag  penetration. 


Volume  Analysis 

(51)  Area  of  Intercept  = AOI 

where, 

X = XBAG  ' ^LINE 
AB 


■/ 


xdy’ 


for  a given  y9  between  yg  and  yA< 
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For  x. 


Mgs 


^G 


= c * 1 


Yl_ 


(Equation  for 
bag  original  shape) 


vf  2 

1 - V dy' 


(52)  BAG  = 


2a 


[y’  V 


a2  - y’2  + a2SIN_1  ( 


,\lyB 


yk 


For  XLINE.-,  using  analytic  geometry. 
Ad 


fXA~  X ' X‘  ' 

XLINEAB  = x.  - I — - y I + 


A \yA  - yB  a, 


:A'*B 
(yA  • V 


y' 


Therefore , 


yB 


(53)  LINE-/ y.x^dy'- 


(*A  - V \ . , (XA  - V 

'A  yA  - yB  WY  2CyA  ' V 


y 


,2 


CHORD 


Under  certain  conditions  (x^  and/or  < 0) , we  have  to  add  " 


ACHORD": 


ACHORD 


, I y* 

= 2 fy  *mG  dy'  = c/a  y'  V a2  ■ y’2  + a2siN  1 Cy'/a)  J 


AOI  = BAG  - LINE  + ACHORD 
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We  now  have  the  terms  necessary  for  the  volume  of  intercept  calculation: 


VOI  = (A0I)WAVG  . 


(54)  VOI  2 W 


AVG  1 2a 


y»  V a2  - y'2  + a2SIN-1^^ 


-1 


mCOSe  - SIN0 
c c 


(COS0c  + mSIN0c)y' 


- By’ 


y’ 


yk 


where 


W 


^FT  " ®BAG^Wb  + RHEAD  WH 


avg  Rpr  ' rbag  + rhead 


SIDE  VIEW  TOP  VIEW 


The  error  introduced  by  adding  in  the  volume  marked  * above  only  occurs  for 
AB  > and  will,  for  this  condition,  be  a very  small  percentage  of  the  total 
volume  which  will  be  compensated  somewhat  by  bag  stretch;  i.e.,  we  may 
approximate  the  missing  volume  * by  (we  assume  LINE  AB  perpendicular  to 
x-axis) : 


W, 


V 


ZhAVG  X 


rxAB 


MISSING 


» 
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where 


Since  we  are  interested  in  this  volume  at  the  worst  condition,  high 

penetrations , let : 


M * 20  inches 
as  - 13  inches 


c*  * 6 inches  (a  bag  26  inches  high  by  1Z  inches  deep) 
3 15  inches 


^MISSING 


~ 13.7  cubic  inches 


V, 


BAG,  ORIG. 


= 4/3tt  a^c  = 4/3ir  (13)  ^(6)  = 4244  cubic  inches 


I ERROR 


13.7 

4244 


0.3%,  a negligible  amount. 


We  therefore  conclude  that  we  may  safely  sidestep  the  computation  of  the 

small  volumes  marked  by  the  *. 
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The  volume  of  the  airbag  for  a given  intercept  AB  is  extremely  difficult  to 
compute  exactly  due  to  the  assymmetry  of  the  volumes  marked  by  + and  - on 
the  previous  page.  In  order  to  facilitate  this  computation  we  make  the 
following  assumption. 

Assumption:  The  bag  volume  for  a given  volume  of  intercept,  VOI,  is  in- 

dependent of  the  torso  inclination  if  the  deformed  periphery  is  the  same; 
i.e. , the  crosshatched  volumes  are  equal  for  constant  periphery  and  constant 
volume  of  intercept  (VOI) . 


The  validity  of  this  assumption  must  be  checked  by  computer  results  versus 
test  results. 

For  subsequent  volume  computations  we  use  the  actual  volume  of  intercept 
for  a link  in  order  to  change  the  as  symmetrical  problem  into  a symmetrical 
one. 


For  the  symmetrical  case  we  assume  a disc  penetrates  the  bag. 
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(58)  VOISYM  - -rra 


(59)  VOISYM  = wa 


. xSYM^  xSYM2 
xSYM  - I — — - — — + xSYM 


xSYM3  xSYM2 


3c 


This  cubic  equation  will  be  solved  in  the  computer  program  by  using  the 
Newton- Raphson  method. 

For  constant  periphery; 

— 

(60)  Ellipse  Perimeter  - 2 (OC)  = 2 BCSYM  + ABSYM 
Rewriting  (60) ; 


2tt  1/  - 2(OC)  = 2 BCSYM  + ABSYM  where, 


ABSYM  = 2y'  @ x*  = c - xSYM  in  — — 


y’  = a 4/1  - 


(c  - xSYM) 


(61)  ABSYM  - 2a  * 1 " 


(c  - xSYM) 


/ \ 


Solving  (60)  for  BCSYM 


(62) 


2tt t/a  y - 2 OC  - 2a 
BCSYM  f 2 


xSYM) 
2 


With  BCSYM  known  we  can  solve  for  RSYM  and  0SYM;  i.e., 
(63)  RSYM  (0SYM)  = BCSYM 
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(64) 

(65) 


RSYM  SIN 


BCSYM 


( 


9SYM  \ 

2 / 


BCSYM 


where , 


= y (x^  - x^,) 2 + (y^  - y^)2  (known) 
We  now  have  2 equations  in  2 unknowns  (63)  § (64) 


ff-fCs  BCSYM  j-yiki 

C66)  w SIN 


(t)  • 


BCSYM 


This  transcendental  equation,  like  Equation  25  will  be  solved  numerically 
on  the  computer  for  RSYM  § 0SYM. 


Now  to  calculate  the  bag  volume. 

x 


(67)  V 


ABCD 


rc 

= tt  / y' 
J c-xS 


2 dx? 

c-xSYM  = Xg 


(68) 


y'  = 


fy'  - y' 

BS  ISl 

x? 

xc 


i*bs 


Substituting  (68)  into  (67) : 


(69)  VABCD= 


yBS  ‘ YC 


A 


- x. 


2y, 


BS 


i*BS 


— 

xc/ 


where  is  the  volume  enclosed  by  the  frustum  ABCD. 
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(70)  VBC  = 2ttRS  RSYM2  (eSYM  - SIN(0SYM)]^  » 


where  V ^ is  the  volume  o£  the  ring  axoimd  the  volume 

RS  = Distance  to  centroid  o£  the  segment  BC  • EC 


(71)  a - TAN 


y*  - y» 

7c 


Hs 


(RSYM 
90  - 

, ___  ___  /eSYM 

s y^  - RSYM  COS  f ■ — — + a 


(73)  RS  = RC  + 4/3 


RSYM  SIN' 


3 / QSYM 


0SYM  - SIN  0SYM 


(74)  VBpiM  = - (c+YctlM1  )^bstituting  69,  70,  and  74  into  75  will 

<b  v A ' 7 yield  the  bag  volume,  V, 


TOTAL* 


(75)  V = V + V + V 
1 j TOTAL  ABCD  BC  BRIM 


For  any  AB  resulting  from  a calculated  volume  of  intercept,  VOI. 


This  completes  the  derivations  necessary  for  programming  the  bag  forces  and 

geometry. 
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APPENDIX  C 


DERIVATION  OF  THE  STEERING  COLUMN  ALGORITHM 


COLUMN  FORCE  CALCULATIONS 


The  total  force  acting  axially  along  the  column  is  given  by: 

(76)  fac  = W “S  (9c  ■ 9t)  + fhead  cos  (ec  • V 

The  total  force  acting  normal  to  the  column  is  given  by: 

(77>  fnc  = W SIN  (9c  - V + fhead  sin  (9c  - V 

The  total  moment  acting  at  point  0 is  given  by: 


(78)  M = F____  • t (t  is  shown  on  page  35  and  calculated  on  page  53 
U Lrihoi 
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GM  Type  Colum 


= reaction  at  shear  capsule 
^BA  ~ reacti°n  at  aft  bushing 
RgP  = reaction  at  forward  bushing 


Solve  for  external  reactions: 


o 

“it 


FW 


£ 


'SC 


NC. 


FW 


The  problem  is  statically  indeterminate;  however,  it  can  be  reduced  to: 


(791  R = F + 
1 J SC  NC 


2 % * FNC  £SC^ 

£fw 
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(80) 


-!(- 


U + F. 


NC 


'FW 


M. 


(81) 


0 


+ F £ 

NC  SC 


We  may  now  solve  for  and  R^. 

aBF  ’ 


3/M0  + FNC  £SC\ 

+ iJ 


+ f £ 

NC  BF 


+ Mo  “ Mo 


+ F. 


NC 


<82>  REA  = 


/M  + F £ 
3/0  NC  SC 
F + —I  

2\  ^ 


^BF  " lS(?  + FNC  (£SC  ‘ £BF^ 


£„„  - £ „ 

BF  BA 


0 = _RBF  ^BF  “ W 


F + 
NC 


3/V 

2\ 


F a 
NC  SC' 


"FW 


^BA  ‘ *SC^ 


F £ — M+M+F  £ 

NC  BA  0 0 NC  SC 
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Solving  for  , 


F + - 
NC  2 


M_  + F 
0 NC  SC 


(32)  - 


JL 


FW 


(aba  " £sc)  + fnc  ^sc  " aba^ 


£bf  " £ba 


Nqte:  For  a pinned  end  at  the  firewall  the  3/2  factor  in  Equations  79,  80, 

82  and  83  is  equal  to  1.0  and  M s 0. 

Solve  for  "t"  the  moment  arm  for  Fp. 

The  methodology  for  this  calculation  is  as  follows: 

1.  Find  equation  for  line  from  x^,y^  perpendicular  to  AB. 

2.  Find  distance  from  this  line  to  x^,y^.  This  distance  is  "t". 

The  equation  for  line  AB  is  given  by  Equation  10,  i.e. , 

B + x8  (m  COS6  - SIN6  ) 

® 'U  Lj 

7 * cose  + m SIN6 

*L  Lj 

where  B s y^  - y^  - m (x?  - x^) . 

The  slope  of  AB  is: 

m cose  - siNe 

j _ La  'Ll 

m ~ cose^  + m SXN0  ~ 

L L 

For  a line  perpendicular  to  this  line, 

^PAB  = " in® 

The  line  from  x^,y^  perpendicular  to  AB  is  given  by 
(y’  ‘ V ’ mPAB  (X'  ’ Xo} 
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Rewriting  the  equation, 


'mpAB  x'  * y'  ~ yQ  + mpy\B  Xo  = 0 
The  distance  t is  given  by 

'mPAB  * yFT  * (-yo  + mPAB  V 
(84)  t = 

V"mb  + 1 

t > 0 for  x^,y^  above  a parallel  line  through  point  0. 
t < 0 for  x^y^  below  a parallel  line  through  point  0. 
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Appendix  D - DRAC  Program  Listing 


> 


- 1 - 


DRAG 


1 £■■■-  0 1 --8 1 


I HOC 

I I OC 
i£oe 

1 3 OC 

1 4 OC 

15  0C 
160C 

1 7 OC 

1 8 OC 
1 90C 
£ 0 OC 
£1  OC 


£4  0C 
£5  OC 


♦♦♦  H PMC  ♦♦♦ 


THIS  PROGRAM  PPEP I CTS  THE  DRIVER  KINEMATICS  im  r CRASH  SITUATION 
IN  WHICH  THE  DRIVER  IS  RESTRAINED  BY  AN  AIRBAG  AND  KNEE  RESTRAINT;  / 
THE  AIRBAG  IS  ATTACHED  TO  A STEERING  CEIL  CNN  AND  WHEEL.  THAT  CELL  APSE'S 
ACCORDING  TO  A PREDETERMINED  FORCE-CRUSH  CHARACTERISTIC. 

THE  KNEE  RESTRAINT  ABSORBS  THE  KINETIC  ENERGY  OF  THE  LOWER  BODY  AND? ' 
LIKE  THE  COLUMN?  CRUSHES  ACCORDING  TO  A PREDETERMINED  FORCE-CRUSH 
CHARACTERISTIC. 

THE  DRIVER  IS  MODELED  BY  THREE  MASSES- A HEAD  MASS?  A TORSO  MASS  AND  . 

A LOWER  BODY  MASS.  THE  DRIVER  IS  CONSTRAINED  TO  HAVE  PLANAR  MOTION 
SO  THAT  THE  PROGRAM  IS  STRICTLY  APPLICABLE  ONLY  TO  FRONTAL  CRASH 
S ITU AT  I DNS . 

EVALUATIONS  OF  AIRBAG?  STEERING  WHEEL?  STEERING  COLUMN?  KNEE- 
RESTRAINT  AND  VEHICLE  PERFORMANCE  CAN  BE  MADE  BY  APPROPRIATE  CHANGES 


£6  OC 
£7  OC 
£8  OC 
£9  OC 
3 0 OC 
31  OC 
3£  OC 
33  OC 


IN  THE  DESIGN  PARAMETERS. 

TYPICAL  DESIGN  PARAMETERS  THAT  CAN  BE  EVALUATED  ARE  BAG  SIZE?  BAG 
SHAPE?  INFLATION  CHARACTERISTICS?  VENT  AREA?  STEERING  COLUMN  AND-- OR 
STEERING  WHEEL  CRUSH  CHARACTERISTICS?  KNEE  RESTRAINT  CRUSH  CHAR- 
ACTERISTICS- STEERING  COLUMN  SUPPORT  STRUCTURE  STIFFNESS*  AS  WELL  AS 
OTHER  SYSTEM  PARAMETERS. 

THIS  PROGRAM  IS  SELF  CONTAINED  IN  THAT  NO  EXTERNAL  FUNCTIONS  OR 
S U B R 0 U TINE S ARE  RE 0 U I R E D . 


34  OC 


35  OC 

AUTHORS  MICHAEL  FITZPATRICK 

36  OC 

F I TZPATR I CK  ENG  I NEER I NG 

37  OC 

W A R S A W ? I N DIANA  4 6 5 8 0 

33  OC 

TEL;  (£19? -£67*4437 

39  OC 

DEC.  I?  1981 

40  OC 

41  OC 

4£  OC 

43  0 

FILENAME  INFILE 

44  0 

REAL  MU  * LSCZ* LFWZ  * LBAZ * LBFZ * KEN 

45  0 

COMMON - OUT -p  NPD  * T (35?  * XO  (6?  35?  * XI  (6*  35?  * X£  (6?  35?  * X3  (6?  35?  ? X4  (6?  35?  ? 

46  0 

&S5  (5  ? 35? 

470 

COMMON  -OUT  I -pH6  ? X6  f 4 * 5 0 ? ? T6  (5  0? 

48  0 

COMMON-' NAME-- INFILE 

49  0 

COMMDN--H I C--TH I C ( 1 75?  * HRGS  1 75?  * P I NT  I 

5 0 0 

I 05  0 

PRINT  ?"  INPUT  FILE  NAME'-’ 

51  0 

INPUT  i INFILE 

5£0 

NPD=  0 

53  0 

H6=  0 

540 

CALL  SOLVE (8? 

55  0 

IF  (NPD. GT. 35?NPD=35 

56  0 

I F (H6 . GT . 5 0?  N6=5 0 

570 

1 1 £ O 

FORMAT  ( I H-? 

580 

1 1 £5 

FORMAT (F7.4* £F6 ; I * £F7 . £ * £FS . 3 ? F3 . I * £F7 . £? 

590 

1 1 3 0 

FORMAT  '.rlX*  7FI I . £? 

54 


BP  MU  It'-'- 01. .'SI 


6 ft  ft  114  0 FORMAT  >:  ) 


6 1 O 

115  0 

FORMAT < IK- 7 ( 47, " 

,•  r 

6£  O 

PRINT  1130 

6 3 1 J 

1155 

FORMAT  C 1 9 1 " 

TIME  CHEST  HP 

CHEST  SI 

HEAD  AP 

6'  4 0 

& HERD  SI "717 

, " C M 8 .)  C b 

" S)  c 

n ' N!  j (,  i n 

65  U 

&:  i.  In  •*  J .i 

66  0 

PRINT  1155 

67  0 

PRINT  1150 

68  O 

DO  1160  K=1,N6 

69  0 

1 1 6 0 

PRINT  1 1 3 0 , T 6 '•  H •' 

< C 3 0’  C J , K > , J = 1 4 4 ) 

7 0 0 

PRINT  1130 

v7 1 0 

117 

0 FORMAT C13n " 

TIME  VEH  b-S 

VEH  VEL 

VEH  DISP 

BODY 

78  0 

& G •*'  S COL  B I SP 

B A b t-'RE.  S S ' ' ■■  ■’  1 X , 

CMS) 

f H •’  S 

CMPH 

73  0 

& CINCHES')' 

CG‘ S)  CINCHES) 

CPS  I b) " 

y 

740 

PRINT  1170 

75  0 

PRINT  1150 

76  0 

DO  1331  K=  1 1 NRB 

770 

1931 

PRINT  1 1 3 0 1 T CK) * 

C 3 0 C 1 4 K ) , J *1,6  > 

78  0 

PRINT  USD 

790 

1363 

FORMAT  Cl 7, " 

TIME  H-P  DISP 

H-P  VEL 

H-P  ACC 

FEM  ED 

S 0 0 

6 RTF  SF AT  PR. 

H-*-P  R : D . 1 3 4 

CMS) 

CINCHES) 

CMPH) 

3 1 O 

S:  C G ‘ S 1 

CLBS)  CLBS) 

CINCHES 

!■  *•'  .*i 

38  0 

PRINT  1333 

33  0 

PRINT  1150 

84  O 

DO  1330  K = 1 , NPD 

85  0 

133  0 

PRINT  1 1 30, T CK) , 

C 3 1 C . J , E ) , J = 1 * 6 ) 

36  O 

PRINT  1130 

37  0 

135  0 

FORMAT  Cl 7, " 

TIME  TORSO  DISP 

TORSO  ANG 

TORSO  VEL 

TORSO 

S3  0 

tv  ALU  TO 0' S D P . D 

TORSO  R. V.  'V 13, " 

CMS) 

CINCHES > 

39  0 

g.:  CDEG') 

C Ii S E C ) C B ■r  S E C ♦ ♦ 3 ) 

CINCHES) 

CMPH)  “) 

9 O i j 

PRINT  I350 

9 1 1 J 

PRINT  1150 

93  0 

DO  1310  K= 1 , NPD 

930 

1 3 1 O 

PRINT  1 130, T CK> , 

•r  X c.  *•  .J  ■»  K 1 * J = I t*  ■l 

94  0 

PRINT  1130 

95  0 

133  0 

FORMAT  Cl 7, " 

TIME  HEAD  DISP 

HEAD  ANG 

HEAD  VEL 

HEAD 

96  O 

& ACC  HEAD  R.D 

HEAD  R.  ANG 13,  " 

CMS) 

CINCHES) 

97  0 

& CDEG) 

C D ••••■  S E C ) C D S E C ♦ ♦ 3 ) 

CINCHES) 

CDEG) " ) 

93  0 

PRINT  1330 

99  0 

PRINT  1150 

1 0 0 0 

DO  138  0 K=I  NPB 

1 0 1 O 

138 

0 PRINT  1 1 30t  T CK)' 

4 C 3 3 C J , K ) <i  J = 1 i*  6 ) 

1 030 

PRINT  1130 

1 030 

140  0 FORMAT  (1^  v 

TIME  COL  A3  FOR 

COL  N FOR 

COL  MOMENT 

COL 

1 04  0 

& RESIST  COL  STROKE  COL  ST  VEL’Vltf, 

CMS 

) CLBS  > 

1 05  0 

& CLBS)' 

C IN-LBS)  CLBS) 

C I NCHE 

3)  C IN-  -SEC 

y " y 

1060  PRINT  1400 

107O  PRINT  115 0 

1080  DO  1460  K= 1 , NPD 

1 0 9 O 1 4 6 O PRINT  11 3 O , T C K ) * C 3 4 C J - K ) - J =1,6  ) 
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DRAG 

!£.-■"  01.  -"81 

1 1 0 0 

PRINT  1120 

1 1 I 0 

148  0 

FORMAT  ,r  I 7 - TINE  ERG  PEN„  ERG 

VOL.  ERG  PRESS. 

i.il  R 

1 1 £ 0 

7 F □ R C E P . F □ R C E " 1 7 - ■ C N 3:  > C I N f:  H E S >'  C C U .IN.  > 

1 I 3 0 

P ,rPS  I G>  CLES>'  CLES>  ">' 

114  0 

PRINT  1480 

1150 

PRINT  1150 

1 1 8 0 

ED  1540  K= 1 - NPD 

117  0 

154  0 

p r i n t 1 1 3 o - T <: k y , c 7 5 < j - k" y , j=i , 5 >' 

1 1 8 0 

PRINT  1 1 £ 0 

1 1 9 0 

183  0 

PR  I NT  1 " ENTER  1 TO  CRLCULRTE  NIC" 

1 3 0 0 

INPUT  1 NRES: 

1 £ 1 0 

I F ■" N RES  . NE.  1 V G □ T □ £ 0 0 0 

1 ££0 

184  0 

PER K = 0 . 

1 £30 

NSTDP=NR 

124  0 

DO  1715  1 = 1 - NS TOP 

1 £5  0 

DO  1718  J= 1 - I 

l£b  0 

L = I + 1 

1 £ 7 0 

SL1M= 0 . 

1 £8  0 

BO  1717  K = I - J 

1 £ 9 0 

L=L- 1 

1 3 1 j 0 

s 1 J M = s u N + HR G s >r L y ♦PINT  1 

1 3 1 0 

1717 

CONTINUE 

.1 3 £ 0 

HELT =TH I C CK > 

1330 

C H E C K = S U N . •••'  D E L.  T 

1 34  0 

I F (PEAK-CHECK > 1 7 1 8 - 1 7 1 8 ? 1 7 1 8 

135  0 

1718 

PFRK=C HECK 

13b  0 

TLOW=  CL-I  ♦PINT1 

137  0 

TNI  GR=  I ♦?’  I NT  I 

138  0 

1718 

CONTINUE 

139  0 

1715 

CONTINUE 

1 4 0 0 

HIC=PERT>*£. 5 

1410 

PR I NT  1 THE  HI C I 5 - H I C 

1 4£  0 

PR  I NT  1 T 1 = " t TLOi-l 

143  0 

PRINT- " T£= " - THIGH 

1435 

£ 0 0 0 

STOP 

1437 

eh  n 

1 44  0C 

145  OC 

THIS:  SUBROUTINE  SETS  UP  THE  B I FFERENT I RL 

EQUATIONS  THAT  DES 

CPI  EE 

146  HI  .J 

THE  H RIVER  KINEMATICS. 

147  0 

S U E P D U TINE  D I F E 0 ( T < V - H Y > 

148  0 

COMMON.--MRNDRT.--ZL  - ZT  - ZH  - RT  - RN  - RH  - RTOPH  - 7£Z 

- YcZ  j i.iJE 

149  0 

B DU  EL  E PRECISION  f C8>' 

1 5 0 0 

DIMENSION  Dt C8> 

1 5 1 0 

C A L. L.  F □ S' C ETH  ■ V - T N E C K > 

1 5£  0 

CALL  DECEL STi GS> 

153  0 

C RL_L  E A G S U E ST  ^ Y - T N F C K - F T H - F S - F T T - G S > 

154  0 

SH=S  I N Cf  <&':•  > 

1550 

S T =S  I N |P Y C7>  >' 

15b  0 

ch=cos  cy  (by y 

157  0 

CT=CDS  CY  C7>"  y 

56 


4 


l!  k'mi:' 


1 if  O 1 .•••‘8 1 


1 58  0 
159  0 
1 6 0 0 
1 6 1 0 
168  0 

163  0 

164  0 
1 6 5 0 
1 66  0 

167  0 

168  0 
1 69  0 

1 7 0 0 
1710 
1780 

173  DC 

174  OH 

175  OG 
1760 

177  0 

178  0 
1790 
1 3 0 0 
1 8 1 0 
188  0 

133  0 

134  0 

185  0 

1 86  0 
1870 
138  0 

1 8 9 0 

1 9 0 0 
1 9 1 0 
1980 
1930 

194  0 

195  0 

1 96  0 

197  0 
1930 
199  0 
8 0 0 0 


DY  ay * -:f czt>rt+2h*rn>' ♦ CCT^DY  <3y-ST^Y  C3>  ♦ Y •' 3 :•  > 

& — Z H ♦ R H ♦ G C H ♦ D Y ,r 8 > - S H ♦ Y •' 8 > ♦ Y C 8 >'  > 1 CZL+ZT+ZH> 

D Y ■: 8 >'  = -I  F T H - Z H ♦ R H ♦ D Y H y ♦CH-ZH+P N ♦PH*  G C T ♦ C H ♦ D Y (' 3 y + C T ♦ S H ♦ Y (. 3 >' 

8 ♦ Y -r  3 y — $ T ♦ F:  H ♦ Y ( 3 :•  Y < 3 > + S T ♦ S H ♦ Ei  Y (.  8 y 1 + Z H ♦ 3 8 . 1 7 ♦ S H ♦ R H ;• 

GZH+RH+EH> 

D Y C 3 y = GFTT-  GZT ♦RT +ZH+EN>  ♦ D Y •'  1 > ♦ G T - Z H ♦ R N ♦ R H ♦ G C T ♦ C H ♦ D Y C £ > + Y C £ • 

8 ♦ Y -r 8 y ♦ C S T ♦ G H - G T ♦ S H 1 +ST ♦ S H ♦ D Y C 8 >'  > +ZT  ♦ ££.  1 7 ♦ R T ♦ S T + Z H ♦ 3 £ „ 1 7 ♦ R N ♦ST >’ 
3 G Z T ♦ E T ♦ P T + Z H ♦ E N ♦ E N y 
DY <4y =-GS 
DY  C5y  =Y  G 1 •- 

DY  *r  6 y xY  t£y 
BY  <77  * Y C3>' 

DY  G8  1 = Y G4 >' 

RETURN 

END 

THIS  SUBROUTINE  READS  IN  THE  INPUT  DATA ? SETS  UP  THE  INPUT  DATA 
FOP  DISPLAY  AND  INITIALIZES  KEY  VARIABLES; 

S U B R 0 U TINE  S E T U P G 7 * Y T 

REAL  MU’  Ml  ’ LSL  ’ LFU’  LBA’  LBF  ’ LSL'Z’  LFWZ ’ LBAZ’  LBFZ » KRN 
C a M M a N s EATFRI  IS N S F , S U N ’ S F U T ’ R E L. S F , S F N C £ - £ 4 > 

C 0 M M 0 N - KNEE R E S T >•  M K R ’ S K R ’ R U T , R E L. K R - K R N ,r£  ? £4  y 
IS  0 M M o N N A M E I N E I L E ’ □ U T F I L E 
E I LENAME  I NF I LE  ’ DUTF I LE 

GOMM0N..-M  AND  A T -- ZL  ’ ZT  ’ ZH » RT  ’ RN  ? RH » RTDF'H  > 7c‘Z  ’ Tc'Z’  UB 
C □ M M 0 N N E G: K N R N ’ FNEG'K  (£•>  84'  ’ DON 
G D M M []  t4  V E H ■-  N V ’ V E H A S G 3 * 8 0 1 
G:  0 M M O N r 5 A S F L 0 r'  N P 6 * A E N f £ h £ 4 i< 

C:  □ M M □ N r is  □ L F □ R --  N P G ’ G 0 L C £’84  y 

is  Q M M a N A A S DAT  A T M D P ’ P G Z ’ G T Z ’ U ’ P N 1 ’ PN 8 P N 3 

COMMON -' BAG  DAT--' VC  1 , VC  8’  AV  » SA,  SC’  71  ’ Y1 

COMMON-- COLDAT  ' SL  I Mu  THETAC ’ MU ’ LSCZ * LFWZ LBAZ » LBFZ ’ WC 

G a M M 0 N WHEEL  I I I H ’ V □ ’ R I M E A D 

COMMON-- PARAU..-RFT’  THETATZ’  THETAHZ 

E 0 M M a N ..-  M I S C .-••  V G 0 L C Q ’ S C 0 L C □ ’ PR 8 

G 0 M M 0 N M M I S C ■-•  F M £ ’ F'A5  ’ G T’FPN’  M G ’ V OLD’  G W » LSC’LFI -I » LBA’  L B F 
DOUBLE  PRECISION  Y G8> 

7=  0. 

£07  0 FORMAT  CV)' 

£ 0 8 0 F D R M AT  a 7 j ■ I HI  T I A L V E L 0 C I T Y i ’ G 1 0 . 3 ■■- 1 7 ■> 

8 • INITI AL  HEAD  ANGLE 3 ■ ’ G 1 0 . 3--'  1 7 » '■  I N I T I AL  TORSO  ANGLE 3 • , 

i:G  1 0 . S i- 
s' 1 0 0 FORMAT f 17’ 


£ 0 1 0 
808  0 

K-  ML EG 

S:  RTOF'H  ■•- 

M TORSO 
1 7 ’ 7 G 1 0 ; ■- 

MAE  AD 

RT 

RN 

RH 

8 03  0 
804  0 

8 1 8 0 

FORMAT  y.  1 7 » ■ 
V,  PM8 

ATMOP 
PN3  ■•-'  1 

PGZ 

7 ’ 7 G 1 0 i 3 .1- 

GTZ 

U 

PH  1 

8 0 5 0 
806  0 

8 1 3 0 

FORMAT  G17’ 

7 1 Z 

VC  1 
Y 1 Z 

VC8 
DCN  ■ •- 

AV 

1 7’  SGI  0;  3y 

SA 

S G 

£07  0 

£138 

FORMAT  a 7’ " 

SL I M 

THETAC 

MU 

LSCZ 

LFWZ 

57 
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DRAG  1 £-'■  Mi  1 -■•'8  1 

£ G B 0 4:  L E A Z LBFZ  W C • 1 ft , 8 G 1 0 ; 3 4 

8 D 9 h i £ 1 3 4 FORMAT  Cl X - ■ i.l H V 0 R I M R h D X 8 Z y 8 Z 

8 1 0 D ft  i .l  B ' ' i x , f,  F i n . 3 y 

8 I 1 U 8 1 35  FORMAT  1 X - " GAS  FLOW  T I ME  “ 1 X - 1 0G I 0 . 34 

8 1 £ O £ 1 3 6 F 0 P M RT  Cl X - 1 ' G R S FLO W - LB •••" S E C ' ‘ 1 X - 1 0 G 1 0 . 3 4 

£ 1 3 0 £ 1 3 Z F □ R M RT  a X - • C 0 L U M N S TRORE  - I N £ H E S *• 1 X , 1 0 G 1 0 . 3 4 

8 1 4 0 £ 1 3 8 FORMAT  Cl X i ■'  £ 0 L IJ M N F 0 R £ E - LB S • •-  1 .X • 1 0 G 1 0 ; 3 4 

£ 150  £ 1 4 0 F 0 R M RT  Cl X - ■ N E £ K R N G L E v 1 X * 1 0 G 1 0 „•  3 4 

£16  0 £ 1 5 0 F □ R M RT  C 1 X « • N E £ K T 0 R © U E.  • 1 X ' 1 O G I 0 . 3 4 ' , 

£ 1 7 O £ 1 6 O F □ R M RTC1 X - 1 • H P T S N E £ K N p T S K R N P T S V EH  HPT S SEAT  H P T S G A S 

£ 1 8 0 ft  H P T S £ 0 L S L ■ S T S L . K R •••  ' 1 X , 8 G 1 0 . 3 4 

£ 1 9 0 £ £ 5 0 F 0 R M RTC1 X - ■ V E H .•  p IJ  LSE  - TIME  M •••" 1 X - 1 0 G 1 0 ; 3 > 

£ £ O O £ £ 6 0 FOR  M R T C 1 X ■,  • VEH=  P IJ  L S E - DEGELERRT 1 0 N • 1 X i 1 O G 1 0 . 3 4 

£ £ 1 0 £ £ 6 5 F □ R M R T c t X - S E R T FRI£TI0N  B I S P L R £ E M EMT  • 1 X , 1 0 G 1 0 . 3 4 

£ £ £ 0 £ £ 67  F0R  M RT  C I X - • S ERT  FRI  £ T 1 0 H FOR  £ E - LB  S • 1 X > 1 0 G 1 O ; 3 > 

£ £ 3 0 £ £ 6 8 F0RMRT  C 1 X - '•  K H E E D I S P L R £ E M E H T • 1 X - 1 0 G 1 0 .•  3 > 

£ £ 4 0 £ £ 69  FO  R M RT  Cl  X - " K H E E F 0 R £ E - L.  B S - 1 X - 1 0 G 1 0 ; 3 >' 

885  0 REHli  CINFILEi.  £ 07  0.-*  V '.4  1 ’ r 1 6 •'  ^ y'  C 7 .-1 

££60  RERD CIMFILEj £070> ZL?  ZT^  ZHu RTj RH. RHj RTGPH 


£87  0 

PERU CIHFILE' £0704 

NF'N  - NK  R i N V i NSF  i NF'G  - NPiS  i SUN  i SKR 

££8  0 

READ CINFILEi £0704 

CGEN Cl < K4 i K*1 1 NPG4 

£89  0 

RERB CINFILEi £0704 

CGEN  C£-  K 4 i K*1 1 NF'G 4 

£3  0 Ij 

PERU CINFILEi £070 4 

CiSOL  Cl , K 4 i K*1 1 NF'G 4 

£31  0 

READ CINFILEi £070 4 

CG0L  C£i  K4  i K * 1 1 NF'G 4 

£3£  0 

READ  CINFILEi £0704 

ATM0P  i F'GZ  i GTZi  Ui  PN1 1 RN£i  PN3 

£38  0 

READ  CINFILEi £070 4 

VG 1 i VC£ i A V - SA i SC i X 1 Z i Y 1 Z i DCN 

£34  0 

READ CINFILEi 807O 4 

SLIM' THE TAG i MU i LSGZ i LFWZi LBAZi LBFZi WG 

£35  0 

READ CINFILEi £0704 

OH  i V0 1 R I MR  AD  ’ XeZ  i V £ Z i WEi 

£36  0 

READ CINFILEi £0704 

CSFN Cl i K 4 i K=1 i NSF 4 

£37  0 

READ CINFILEi £0704 

CSFN C £ n K 4 i K=1 1 NSF 4 

£38  0 

READ  CINFILEi  £070.' 

CFNEGK  Cl  ' K 4 i K*1  i NF'N 4 

£39  0 

READ  C.INFILE4  £070 4 

CFNECK  C£,  K 4 i K*1  i NF'N 4 

£4  0 0 

READ CINFILEi 8D704 

CVEHGS  '•  1 1 K 4 i K*  1 1 NV4 

£41  0 

READ CINFILEi £0704 

CVEHGS C£ , K4 iK*li NV4 

£4£0 

READ  CINFILEi £0704 

CKRN C I , K 4 i K* 1 1 NKR4 

£43  0 

RE AD CINFILEi £07 04 

CKRN  C£i  K4  iK*li  N PF'4 

£44  0 

MG  0=1 

£45  0 

PRINT  £490 

£46  0 

PRINT  i1' INPUT  VAL 

LIES  — INPUT  UNITSC  MSEGi  MF'Hi  DEGREE 

£47  0 

& INCHES*  LBS i FT 

-LBS'  G-S4‘‘ 

£48  0 

GO  T0  £5£  0 

£49  0 

£480 

FORMAT ClXi 1 0G1 0. 3 

4 

£5  0 0 

£49  0 

FORMAT C1H-4 

£510 

£5  0 0 

P R I NT  i I N I T I A L V A L U E S - L 0 N V E R T E D U N I T S C S EC i FT S E C i 

£5£  0 

ft  RADI ANSi  FTi  LB 

Si  FT-LBS i FT.--SEG^^£4  ■ 

£53  0 

£5£  0 

PRINT  £080i Y C 4 4 ^ r 

C6  ■'  i V C74 

£54  0 

PRINT  £ 1 O 0 1 Z L i Z T i 

ZHi RT i RNi RHi RTGPH 

£ 5 5 O P PINT  £ 1 6 0 - N P H ^ H K R - H V ? H S F - H P G ' H P £ - S IJ  H ^ S K R 

£56  O PR  I NT  £ 1 35  < CGEH  C I » K>  ? K*  1 , NF'G> 

£ 5 7 0 P PINT  £ 1 3 6 - C G E N C £ - K > , K*  1 * H G > 
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Ei+'RG 

1 £+01 

+81 

£58  0 

PRINT 

£ I 3 7 1 CuDL  '•  1 - E a K—l » 

MPCT 

CL  y U 

PRINT 

£ 1 38  ’ CL  DL  ,r£  i K .>  9 K—  1 » 

NPC> 

£6  0 0 

PRINT 

££65  CSFN  >U  K>  - K=t  ' 

NSF  > 

£6 1 0 

PRINT 

£ £ 67  9 C x F M C £ i K ) i E = I ? 

N 3 F 

£6£  0 

PRINT 

£ I 4 Li  i C F N E ij  E C 1 ’ E j i E = 

I ' NPN  > 

£63  0 

PRINT 

£15  0 h CP  NECK  C£k  k>  « E'  = 

1 * NPN  :• 

£64  0 

PRINT 

££50’  CVEHGS  '"1  - K>  , K = 

1 ' NVX 

£65  0 

PRINT 

c c U * '.r  V fc  H l ? y c.'.  ? K -1  * (■•'  — 

1 ' NV  ? 

£ 6 6 0 

PRINT 

££68 1 CKRN Cl , KJ  , K=1 1 

r j e r > 

£67  0 

PRINT 

c'  c t' ? v K f?  N '•  c * K 1 * K = I f 

NEP> 

£68  0 

PRINT 

£ 1 £ 0 h R T M E! P * p G Z i*  GTZ* 

ij , PN 1 ' PN£ ' PN3 

£69  0 

PRINT 

£ 1 3 0 !•  VC  1 9 'vG£  i RV  9 SR  ■> 

U ’ 7 I Z ' T 1 Z ’ Ii  L N 

£7  0 0 

PRINT 

£ 1 3£  9 SL I Nn THETRC - NU 

' LSCZ ' L F 1.0 Z ' LBRZ'  LBFZ?  1.0 1 

£710 

PRINT 

£ 134'  i.O H i t □ h R I MRRD ' 7 

£Z  9 i £ Z 9 1.0 Ei 

£ 7 £ 0 

GO  TD 

C£7 1 0 ’ 3 I 3 0 O ' MGD 

£73  0 

£710  PRINT 

£49  0 

£74  0 

V £ > = 

ij „ 

£75  0 

Y C3’>  = 

0 . 

£76  0 

7 ,r4  :'  = 

y‘  C 4 ■'  +1  „ 4666 6 €< 7 

— i — i 

? ,r  cr  , 

lL  i*  i*  IJ  • 

1 '•  .J  ' — 

IJ  0 

£78  0 

'T  G 6 = 

T C 6 > ♦. 01745 3 £ 9 

7 9 1 j 

TRET  RH7 

8 0 0 

‘f  C7  y =7  c 

8 1 Li 

THETRTZ 

£83  0 
£84  0 
£85  0 
£86  0 


£9  0 0 
£91  0 


0 

0 


£ 9 3 £ 


£94  0 
£95  0 
£ 9 6 U 
£97  0 
£98  0 
£99  0 
3 0 0 0 
3 0 1 0 
3 0£  0 
3 0 3 0 
304  0 


y 'r  I > = y ,r4  • 

ZL=ZL  ••"3£ . 17 
ZT=ZT.--3£.  17 
ZH=ZH.--3£.  17 
RT =RT.--'  1 £ . 

PN=PN--  1£. 

RH=E‘H.-'  1 £ . 

RTDPH=RTDPH  - 1 £ . 

3 i.  1 M = 3 ! J N ♦ l £ . 

SKR*SKR^1£. 

THETRC  :*THETRC . 0 1 745389 

RDD=  8 G ♦ 3 9 P T C 1 : - R I M P fl D * * £ 7 R ♦ ♦ ,3 > 

7 1 * 7 I Z + R D El  ♦ G □ 3 CTHETftG  > 
r 1 =■=  7 1 Z + R D D * 3 1 H C T H E T R G > 

VGDLLGD* 0 . 

SGDLGD* 0 . 

PR8=  '■"£ . CPN  1 + 1 . > ♦♦CPN1  CRN  1 - 1 . > 

F M £ = 'v1  C £ ♦ S 9 R T CPR8++  f£ . .--PM  1 -PP8++  c CPN 1 + 1 . > - PM  1 ? 

PR=PGZ+RT  ME3P 

PR5=PR 

GT =GTZ 

FPM=PM£ 

MC:=iilG  386 . 

V □ L Z = 4 . 3 . ♦3.  1 4 1 5 9 ♦ 3 R ♦ ♦ £ ♦ 3 C 
V □LD=‘v  DLZ 
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DRAG 


1 0 1 "8 1 


3 05  0 


3 0 7 0 
3 1 j 0 fi 
3 09  0 
3 1 0 0 
3 1 1 0 
3 1 3: 0 
3 .1  3 0 

3 1 4 0 

3 1 5 0 
3 1 G i J 
3 1 7 0 

3 1 8 0 

3 1 9 0 
30  ij  0 
3010 


30:3  0 

304  0 

305  0C 
308  OC 


oc 

0 


309  0 


3 3 0 ij 
3310 
330  0 
3 3 3 1 J 
334  0 
3 3 5 0 


330  0 
337  OG 
3 38  OG 
339  OG 


34  0 0 
3410 
340  0 

34  3 0 

344  0 

345  0 
3 4 0 0 

347  0 

348  0 
3 4 9 0 

35  0 0 
351  0 
350  0 

353  OG 

354  OG 


G I n = ( p A ♦ V n L Z >■  U ♦GTZ  ? 

LSC=LSCZ 

LFW=LFWZ 

LEA=LBAZ 

LEF=LEFZ 

DO  3010  J^ 1 , 84 

AG  A 1 1 ’ J 1 =0EA  ’■  I * J 1 .■■■  1 0 li  0 >- 

FAEGK  (.  1 j J>  *FNEGK  ( 1 - J?  ♦.  01 745309 

3 FA  (1 , J V =3 FA  f 1 i J>  .-10. 

3 ij  1 0 GPA  •"  1 , ,J  :<  =KRA  •:  1 J 1 0 . 

DO  3050  J=Ii30 

'v1  E H 6 S 'r  1 1 J V = V E H A S 1 - J ? 1 0 0 0 . 

3 1 j 5 0 V E A 0 S f 0 , J j = V E H G S «r  £ - J 1 ♦ 3 0 . 1 7 
SFUT = 0 . 

RUT  = 0 . 

AGL1=0 

GD  TU  £500 
3130  GD AT I AUE 
RETURA 
EAD 

THIS  SUE ROUT  I AE  IS  A GEAERAL I ZED  TABLE  LOOKUP  A AD  I ATERPDLAT I DA 
RDUTIAE  WHICH  IS  GALLED  BY  OTHER  ROUT  I AES . 

3 U B P D U T I r j E L D D G IJ P 'rA,  F U A - A P T S - E ? 

D I ME  AS  I DA  F IJ  A 00i  3 0? 

DD  3.19  0 J = 1 1 ARTS' 

3 1 9 0 I F >'F  IJ  A *r  1 1 J ? . G T . A ? GOTO  3 0 0 0 
30  0 0 IF  CJ. EO. 1 ? J=0 
G=J- 1 

Ei31  (A-FUA  ,r  I i G?  ? ♦ (.rFUA  f0i  J :•  — FUA  '!'0i  G?  ?.••■'  ,rFUA  ■!  1 1 J?  —FUN  (.  1 1 K?  ? +FIJH  •■£ » G ? 

RETURA 

EAD 

THIS  SUBRDUTIAE  GALGULATES  THE  AEGG  TORQUE  AS  A FUACTIDA  OF  THE 
A EGG  AAGLE . 

S U E R D U T I A E F D R G E T H >.r  Y i T A E C K ? 

GDAADA.- AEGK.-'  WF'Ai  F NEE  G r£,  £4  ' , DGA 
DOUBLE  PRECIS  IDA  Yt8? 

TNEGG*0i 

TBAARar-DGA*  •' Y (£?  - Y C3?  ? 
vREL=Y  f£? -Y  C3? 

TREL*Y  'IE?  -Y  C7? 

IF  -:TREL.  GT.  0.  0.  AND.  VREL.  LT.  0.  ? GD  TO  I 0 
I F '"TREE . L.T . 0.  0 . A AD . VREL . GT . 0 . ? GO  TO  1 0 
GALL  L DOG  UP  'ITREL  i FAEGG  i AF'A , TAEGG? 

1 1 j T A E C = T A E G G + T D A M P 
RETURA 
EAD 

T H I S S U E R 0 IJ T I A E D B T A I A S THE  C R A S H p IJ L S E G - S A S A F U A C T I DA  OF'  TIME. 


60 


DR  AC 


1 £-■-  h 1 .••-8 1 


355U 

3560 

3570 

358  0 

359  0 
36  0 0C 
36 1 OC 

368  OC 

36  3 OC 

364  OC 

365  OC 

366  OC 

367  0 
363  0 

369  0 

37  0 0 
371  0 
3780 
37  3 0 

374  0 

375  0 

376  0 


379  L 
38  0 
38 1 L 
388  L 
333  i 
384 1: 


387  0 

388  0 

389  0 
39  0 0 
39 1 0 
398  OC 

393  0 

394  OC 

395  0 

396  OC 

397  OC 

398  0 
3990 
4 0 0 0 
4 0 1 0 
4 08  0 
4 03  0 
4 040 


S U E R □ d T I N E D E C EL  ‘T , G S 
C □ M M n N v E H N V , V E H G S '■  £ - 3 0 
C ALL  L □ □ K U RCT,  V E H G S - N V - G S >' 

RETURN 

END 

THIS  SUBROUTINE  COMPRISES  THE  MAJOR  PART  LIE  THE  DRAG  PROGRAM.  IT 
EVALUATES  THE  BAG  SHAPE  AS  A FUNCTION  OF  BAG  PENETRATION  AND  TORSO 
ANGLE - CALCULATES  THE  FORCES  THE  BAG  APPLIES  TO  THE  DRIVER-  CALC- 
ULATES THE  BAG  VOLUME  AND  PRESSURE-  DETERMINES  THE  GAS  GENERATOR 
FLOW  CHARACTERISTICS  AND  CALCULATES  THE  STEERING  COLUMN 
FORCES-  MOMENTS  AND  STROKE ; 

S U B P OUT  I ME  BA G S U B ( 7 - V - TNEGK  - FTH  - F 7 , F T T - G S >' 

REAL  M P U - M U - L S L - L. F W - L BA-  L B F - L F W 7 - L 3 C 7 - LBAZ  ? LBF  Z j M 0 - M L 
REAL  KPN 

G 0 M M 0 N S EATFRI C N S F - S U N - S F U T ? PEL S F - S F M C 8 - 8 4 
C 0 M M 0 N M A N DAT Z L - ZT  - ZH-  RT  - RN-  RH-  RTOPH-  78Z-  Yc'Z?  WB 
C 0 M M 0 N K N E E R E S T N K R - S K R - R U T -PEL K R - K R N t £ , 8 4 
0 0 M M 0 N G A S F L 0 - N P G - G E N ( £ -84' 

C 0 M M 0 N ••••  c 0 L F 0 P ■•■■■  N P C - C D L.  C £ - 8 4 > 

C 0 M M o N G A s D A T A 7 M 0 P - P G Z - G T Z - U - F'N  t - F‘ N £ - P M 3 
C 0 M M 0 N B A G D A T --' V C 1 - VC 8-  AV-  SA-  SC  - 71  - VI 

C 0 M M D N , - COL.  DAT  SLIM-THETA  C - Dt  U - L S C 7-  LF  W Z- LBAZ- LBF Z-  W C 
C 0 M M 0 N W HEEL Dl  H - V 0 - R I M R A D 
C 0 M M 0 N •••••  M I S C V C 0 L C 0 - S C □ L C D-PR 8 

COMMON.-' MM  I SC.--FM8  - PA5-  GT  - FF'N-  MC  - VDLO  - GW  - LSC  - LFW  - LB  A - LBF 
COMMON-- PARAM-  RFT-  THETATZ-  THETAHZ 
C: 0 M M 0 N TIME S TEP-  7 S T 0 F' 

C: 0 M M 0 N M PAPA r'l F A C 0 L - S C 0 L C - P G 1 - F K N EE  - SF-  F N C 0 L.  -MG-  F R C □ L - V C 0 L_ C 
is Q M M 0 N M f'1  PAPA M --  B P - V 0 L - F F T - F F' 

DOUBLE  PREG I S I ON  V C8> - B - A - D - E - A 1 - B 1 - C 1 - 78  - V£ 

5 FORMAT  OH-) 

DIMENSION  DV  .r8>' 

1.1  B ACT =1.1  B 
WHACT =WH 

thet  at =v  <:?y 

THET AH=V  'r6> 

CHECK  TO  SEE  IF  DRIVER  SUBMARINING; 

IF  CAES  (THETAC-THETAT? ; GT. 1 ; 4>  GO  TO  500 
CALCULATE  THE  SLOPE  OF  THE  DRIVER  TORSO; 

□ SLDPE=T AN  y 3 . 1 4 1 59-- £ . +THETAT> 

CALCULATE  THE  NEW  H-F'OINT  COORDINATES  AND  THE  7-COORDINATE  OF  THE 
POINT  WHERE  THE  RIM  INTERSECTS  THE  BAG; 

7£*7£Z-  y'i  .:5>  -V  C8T > ♦ 1 8 . 

V 8 = V 8 Z 

7 IS  * - S IS  ♦ S 0 R T n ; - R I M R A D ♦ ♦ £ S A ♦ ♦ £ > 

B = V 8 - Y 1 - 0 S L 0 P E ♦ C 7 £ - 7 1 > 

A=SA**8+  AS  03  C THET  AO  +OSLOPE^S  I N ,r  THE  TACT  > ♦♦£ 

Z, + S IS ♦♦£♦  -r S IN  CTHETA C >' - 0 S L 0 P E ♦ G OS  '"THET A C T > ♦♦£ 

D=8;  ♦B^SC^^S^  (0  SLOPE  ♦COS  THET  AO  -SIN  •:  THET  AO-  T 
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DRRG 


4 0*50 
4 06  06 
4 07O 
4 08  06 
4 09  0 
4 1 0O6 
4 1 1 0 
4 1 £ 0 
4 1 3 0 

414  0 

415  0 
4 1 6 0 
417  0 
4 1 8 0 
4 1 9 0 
48  0 O 
481  06 
438  0 
483  06 
4840 
4 £ 5 0 
486  0 
4870 

438  O 

439  0 
4 3 0 O 
431  O 
438  0 
4 330 
4 3 4 O 
435  0 
4 36  0 
4 37  0 
438  0 
4 39  0 
44  0 0 
44  1 0 
443  0 

443  0 

444  0 

445  06 

446  0 

447  06 

448  0 

44  9 0 

45  0 06 
451  O 
458  0 
4530 
454  0 


— M — 


1 £■■••'  0 1 •••"8 1 


E*B**£*  65G**8>  -SR**£*  63G**8>  ♦ (60S  '-THETRG  :■  +0SLOPE>SIN  '"THETRGX  X ♦ ♦£ 
TEST  F0P  SIGH  06  B I S6R I H I NRTE 

6 IF  ,r 6) ♦ ♦ c! — 4 „ ♦ R ♦ F .i  y o j 7 !•  7 

RERL  DISTINCT  R00TS  '"DEFINITE  T0RS0  RND  ERG  GONTRGTX 

7 DISG* CD*^£-4. ♦ R ♦ E X ♦ ♦ ; 5 

ERG  INTERCEPT  POINTS-  7R-7B  RND  TR * YE 
7 R*  C-D- D I S 6 X C £ ; ♦RX 

I F 656  „ LE . RES  67 R>  y 7R=RBS  67  R X 67  R*  '"  36- . 0 0 1 X 
7E*  '"-D+DISC  6 C3.  ♦ RX 

I F '"SC . LE;  RES  67 E > X 7B*RES  67 B > 67 E*  CSC:  - . 0 0 1 X 

Y R*  <E+ 7 R*  CO  5 L 0 F* E ♦ 6 0 3 -:TRFTR 6 X - S IN  ■:  THETRG  X :•  X 6 C 6 QS  CTHETR 6 :•  + 0 3 L. □ P E ♦ 
!i.:S IN  C THETRG XX 

I F CS R.  LE . RES  CYR  > YR*RBS  67 R1  6YR-*-  65 R- . 0 0 1 

Y E * CB+ 7 E ■ ♦ C 0 3 LOPE*  6 Q 7 CTHETR 6 >' -SIN  CTHETR 6 X X X 6 C 6 0 3 CTHETR 6 X + 0 S L 0 P E ♦ 

;7 SIN  CTHET RGX  X 

IF  CSR „ LE;  RB 3 CYB :•  X YB*RBS  CYBX  6 YE:*-  CSR- ; 0 0 1 X 
RESTED I STRN6E  FROM  P0INT  R TO  POINT  E. 

R B 3 T 3=  5 Q R T C C 7 R - 7 E X ♦ ♦ £ + C Y R - Y B X * * 3 X 

7 RND  Y CODED;  0F  H-P0INT  IN  7F‘P  I ME  - YPRIME  600RD.  SYSTEM 
7H*  '"  CY8-Y 1 X * 5 1 N CTHETR6>'  + '758-7 1 X ♦60S  CTHET RGX  X 

Y H*  C C Y £ - Y 1 V - C 7 £ - 7 1 7 ♦TRN  CTHETR6  X X ♦ C + 6 0 3 C THETRG:  X X 
IF  CTHETRT-THETR6':'  8?  8-  9 

8 YP*YB 
YN*YR 
7F'*7B 
7N*7R 

GO  TO  ID 

9 7 F'  * Y R 
YH*YB 
7P*7R 
7M*7  E 

1 0 k B R 6 * N F1  T C C 7 N — 7 H 1 ♦ ♦ y + 1 7 N — t H .•  ♦ ♦ y 
I F CREST +RBRG- 1 £ . ♦RN>  11-  1 1 - 18: 

1 1 7FT*  C7R+7EX  •"'£. 

Y FT=  C Y R + Y B 8 ■ 

GO  TO  13 

18  7NE6K*7N-  CRN*  18.  -RERGX ♦SIN  CTHETRT-THET RGX 
YNE6K  *YN+ CRN* 1 8. -RERG> ♦6GS CTHETRT-THET  RGX 
7 FT* C 7N+ 7N E 6 KT 68; 

YFT* CYN+YNEGTO  63. 

RFT*D I STRN6E  FROM  H-P0INT  TO  POINT  OF  FORCE  RPPLI6RTI0N  ON  TORSO 
13  RFT* S Q R T C C 7 H - 7 F T X ♦ ♦ 8 + C Y H - Y F T X ♦ ♦ 8 X 
SLOPE  OF  LINE  PERPEND I 6ULRR  TO  RE 

P SLOPE*- CC  DS  CTHET RGX +0SL0PE*S I N CTHETR6XX8 COSLOPE^GOS CTHET RGX 
3-SIN  CTHETR 6 XX 

T* MOMENT  RRM  OF  TORSO  FORGES 
70*- SC: 

T * C - 7 F T ♦ P S LOPE+ Y FT- Y 0 + P S L 0 P E ♦ 7 0 X 6 S 0 PTC  P SLOP £♦♦£+!; X 
R 1 * 1 . 6 S'  6 ♦ ♦ 8 + P S L 0 P E ♦ ♦ £ 6 S R ♦ ♦ 8' 

E 1 * 8 ; ♦ P S LOPE 6 S R ♦ ♦ £ ♦ C Y F T - P S L D P E ♦ 7 F T X 
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12.-U1.-S1 


- 1U 


Ei  Phl 


455  0 

456  0C 

457  0R 

458  0 

459  0 
46  00 

461  0 

462  0 

463  0 
4640C 

4 6 5 U 

466  00 

467  0 

4 6 8 0 
469  00 

47  00 
471  00 
4720 

473  00 

474  00 
4750 G 
476  0 6 
47700 
4780 
4790 
43  0 0 

431  0 

432  0 

433  0 

434  00 

435  0 

436  OS 

487  00 

488  00 
439  0 
49  0 0 

491  0 

492  0 

493  0 

494  0 
4950 

496  0 

497  0 
4980 
499  0 

5 0 0 0 
5 0 1 0 
502  0 
5 0 3 0 
5 04  00 


16 


17 


18 


21 


0 1 = C Y F T - P S L 0 F’  E ♦ 7 F T Y ♦ ♦2.--'SR^2-  1 . 

PD  I NTS  G RNE'  I HPE  THE  PD  I NTS:  WHERE  THE  LINE  DF  RET  I DM  DF  THE  FOPEE 
l.il  DDL  Ei  INTERSECT  THE  ERG. 

7 G*  C-B1  + S Q P T C E 1 * * 2 - 4 ; ♦ R I ♦ 0 1 >'  Y C 2 ; ♦ R 1 >' 

7 I * C-B 1 -SQPT  CB 1 ♦♦2-4 ; #R  1 ♦ 0 1 >'  Y .--  C 2 . ♦ R I Y 

Y G * P S L D P E ♦ C 7 G - 7 F T Y + V F T 
V I = P S LDPE*  C 7 1 - 7 F T Y + V F T 

Y P I M = R S L D P E ♦ C 7 0 - 7 F T > + Y F T 

G I * vwFT  C C 7 I — Y:  G ■'  ♦♦c!+  1 'T  I — t G '•  ♦ ♦EY 
0RL0HLRTE  THE  ERG  F'ENETPRT  IDN. 

E F' 1 L'-1  F' T ( C 7 FT- 7 G y ♦♦  2 + C Y F T—  Y G Y ♦ ♦ 2 Y 
BETEPMINE  THE  MIDPOINT  DF  LINE  GI. 

7 M G 1 — C 7 G + 7 1 .•••'  2 . 

YNG I — ,rYG  + Y I Y •••■'2 . 

0RL0DLRTE  MRJDP  R7IS  LENGTH  DF  ELLIPSE  PEPF'ENE I GULRP  TD  TDPSD 
ZP* S 0 P T C S R ♦ ♦ 2 ♦ Cl—  7 M G I ♦♦  2 ••••' S G- ♦♦2 •'  — Y N G I ♦♦  2 >' 

0RL0DLRTE  PEP  I NETEP  DF  ELLIPSE  PEPF'ENEi  I GULftP  TD  TDPSD 
P E P B B =«=  2 .•  ♦ 3 ; 1 4 1 5 9 ♦ S 9.  P T C -r  C G I 2 .•  > ♦ ♦ 2 ■ + 2 P ♦ * 2 > 2 . > 

SOLVE  PDF'  THE  RNGLE  CEETR > THRT  THE  F REP  10  TENSION  FDP0E  0DNPDNENT 
MRKES  WITH  PESPE0T  TD  R LINE  NDPNRL  TD  THE  TDPSD;  FIRST  SOLVE  FDP 
F'H  1 1 U S I N G THE  N E l.l TON  PRPH S 0 N M ETH 0 B DP  S D L. V I N G TPR N 0 E NBENTRL  E 9. N S , 
LET  PHID*RN  ESTINRTE  DF  THE  PDDT  PHI  RNB  EPSLDN*THE  DESIRED  RCC- 
UPR0Y  DP  THE  PDDT. 

PH ID-3. 14 
EF'SLDN-;  0 0 001 
PH  I -PH  ID 

I F CRES  CYP I M :<  -P I MPRB > 1 7 j 1 6 - 1 6 
P I M- 0 . 

GDTD13 

■PIN-  IS  THE  0HQPB  LENGTH  DF  THE  PIN  RT  POINT  I. 

P I M - S 9 P T C P I N P R D ♦ ♦ 2 - Y P I M ♦ ♦ £ 1 ♦ 2 . 

PDF'  ERG  F’ENETPRT  I DNS  LESS  THRN  DNE-HRLF  THE  SHE  ST  THI0KNESS-  THE 
EDDY  WIDTH  IN  EON  TRET  WITH  THE  ERG  WILL  NOT  E.7EEED  THE  LENGTH  DF 
THE  EDDY  IN  ED NT RET  WITH  THE  ERG. 

0 H E S T T - l.il  E R E T 2 . 5 

IF  CEP.  LT.  0HESTT.-2.  . RNE;  WERGT.  GT.  REST  - WE-REST 
EBP  - 39 PT  C 030-7 1 :•  ♦♦2+  CYP  I M- Y I >'  ♦♦2Y 
S E P - S 9 P T C ; 2 5 ♦ P I M ♦ ♦ 2 + E E P ♦ ♦ 2 Y ♦ 2 . 

FPHI-CPEPEE-  1.1  E - S E P > ♦ s I N C P H 1 0 2 ; Y PH  I -CGI-  E P - B E P Y C C:  D S CRTR  N C C W E-PI  M Y 
i:  / C 2 ; ♦ C G I — El  F'  — El  El  P Y Y Y Y Y 

EFF'H  I - CPEPEB-WE-SEPY  ♦EDS  CPH I .-"2 . Y C2 . ♦F’H  I Y - C CPEPBB-WB-SBPY  ♦S I N CF'H  I .--2 
Y F‘  H I ♦ ♦ c!  Y 
DEL—FPHI.--EFPHI 


PHI-F'HI+DEL 

I F CF’H  I ; GT ; 2 ; ♦S ; 1 4 1 593 ; DP ; F'H  I . LT . 0 . Y GD  TD  52  0 
IF  CRES:  CBELY  . LE;  EF'SLDNY  GD  TO  22 
GD  TD  21 

22  RLPHR-RTRN  C CWB-PI N Y C 2 ; ♦ C G I-EF’-D E P Y Y Y 


EETR-PHI.--2;  +RLPHR 

SDLVE  FDP  THE  RRE  SIN  DF  YF'.--'SR  RNE  YN  - SR . 


63 


1 c' ■-Hi  ■■-SI 


DR  fill: 


5 05  0 


5 06  0 
5 07  00 
5 08  06 
5 0 9 0 


5 1 0 0 
5 110 
5 1 £ 0 
5 1 3 0 
5 1 4 0 
5 .1 5 0 

5 1 6 0 

517  0 
51800 


519  0 


5£  0 00 

5£1 

0 

5££ 

0 

cr  ~i  “i 

j cl 

00: 

5£4  00: 

5c,c5 

0i_: 

526 

0 

5? 

5l 

50 

5' 

5:' 


8 0 
9 0 
0 0 
1 0 
£0 


£5  HSYPSfi*fiT fiW  •:  ( YFvSfi  :•  (ISQPT  C 1 . - ( Y P S fi  > y y 
fi S YN  Sfi*fi  TfiM  1 (YH/'Sfi)  .•••'  CSQPT  Cl ; — (.  YN  ••'SA.'1  ♦♦£>  > > 

BAD-  TLINE  AND  flCHDRD  fiRE  I NTERMED I fiTE  VALUES  REQUIRED  FDR  THE 
fiREfi  DF  I NTEREEF'T  CALCULATION; 

B fi  b * S 0:  .••■•  ; ♦ 3 fi  :•  ♦ C C V P ♦ S Q R T C S fi  ♦ ♦ £ - V P ♦ ♦ £ > + S fi  ♦■*•£♦  fi  S Y P S fi  > - 

& C Y H * 3 Q E T r S:  fi  ♦ ♦ £ - Y N ♦ ♦ £ :•  + S A ♦ ♦ £ ♦ fi  S Y N S fi  >'  > 

TE I ME*  CXA-  CXfi— XB.:*  ( Yfi- YB>  ♦Yfi.:'  ♦YP+  CXA— X By  •••■'  C£;  ♦ v‘ fi — YB  -1  > ♦ YP ♦♦£— 

& ( Yfi-  ( Xfi-X B:  (Yfi- YB.:*  ♦Yfi.:*  ♦YN—  (.Xfi— XB-’1  (£;  ♦ (Yfi- YE:1  .’|  ♦ Y !'!♦♦£ 

fiEHDR D^ 0 „ 

I F C X P;LT.  0 . > fi  ‘3  H □ R B * S C •••"  S fi  ♦ C 3 ; 1 4 1 5 9 3 ♦ S fi  ♦ ♦ £ £ . - C Y P ♦ 3 Q R T C S fi  ♦ ♦ £ - Y p ♦ ♦ £ y 
R- + 5 fi  ■*•♦£♦  fi  S Y P S fi  > :< 

IF  C X N„  LT . 0 ; > fi  0:  H 0 RB*fiC'H  □ R D + S'  0:  .•••'  S fi  ♦ C Y N ♦ S Q R T C S fi  ♦ ♦ £ - Y N ♦♦£  + S fi  ♦ ♦ £ ♦ fi  S:  Y N S fi 
8 + 3 ; 141  5 9 ■ ♦ S fi  ♦ ♦ £ £ „ > 

SOLVE  FDR  THE  fiREfi  DF  I NTERC EFT . 

ADI =Efib-TL I NE+fiGHORB 
■ V □ I - * V □ L U M E D F Bfi b I NTEREEF'T . 

i i ! fi V b ^ 1 >.  kF T — PEfib  1 ♦WB+RHE AB+WH,1  ■■  * P1 1 — F — R B h b R H E H D •* 

V D I = i.il  fi  V b ♦fiDI 

THE  FDLLDi.il Mb  ROUTINE  USES  THE  NEWTDN-RfiPHSDN  METHOD  TD  SOLVE 
fi  EUR  IE  EQUATION  FDR  THE  Efib  PENETRATION  THfiTMDULD  EXIST  FDR 
THE  3 Y M M E T R I E fi  L C fi  S E WITH  fi  b I V E M V D I . 

RDDTb*0.-  8+ BP 

I F ( B P . b E . S E +AB S'  ( XI U ■'  ■'  ROOTb* S' U 
RDDT*ROOTb 

7 3 FROOT* 3 . 1 4 1 5 9 ■ ♦ S fi ♦♦£♦(- R DO ' T * ♦ 3 C 3 . ♦ S E ♦ ♦£>  +RDDT ♦ ♦ £ S 0 > - V 0 1 
D F R □ 0 T *3.  141 5 9 * S A ♦•*•£■*■  <— R D DT ♦ ♦ £ S E ■*■♦£+£ ; ♦ROOT S E > 

DELRT  ^-FPDDT  DFRDDT 
ROOT =RDDT+DELRT 


IF  CPODT. bT. £. ♦SE. DR. 

534  0 

IF  CfiES ( DELFT ' . LE. EP 

5 9 5 0 

bD  TD  73 

5 3 6 0 

74  PEE=£ . ♦3  ■ 141 59+SQRT 

537  0 

BPSYM*RDOT 

538  0 

fiBSYM*£ i ♦Sfi+S QRT C 1 . - 

539  O 

YE=R I MR AD 

54  O O 

DC*SQRT ( CXC-XOT ♦♦£+ ( 

54.1  O 

B G: S Y M E = ( PER— £ . ♦OC -AB 

54  £0 

XES=SE-EPSYM 

543  0 

YBS-fiBSYM  -'£. 

544  OC 

545  OC 

546  oc 
547U 

548  0 

549  0 
55  O 00: 
551  OC 
5 5 £00: 

cr  cr 
! . I . I I 


T . 0 „ > bD  TD  518 
DO  TO  74 

+ S E ♦ ♦ £ >"  8 • :• 


BEFORE  C ALGOL AT  I Nb  THE  LENbTH  DF  THE  LINE  EE  FDR  THE  SYMMETRICAL 
EASE  CHEEK  TD  SEE  THAT  THE  HEAD  AND  CHEST  HAVE  NOT  EDTTDMED 
OUT  ON  THE  WHEEL  RIM.  IF  THEY  HAVE-.  STOP  THE  RUN. 

I F CXfi .LT. XE . AND . Yfi . LT. 0. > bD  TD  5 1 0 
I F C IX fi . L T . X E . fi N Ei . Y A . b T • U . .:•  bO  TD  518 
B i..:  s Y M x. = x.  i.;jR'  T C C X B s ~ X L .:•  ♦ ♦ 8'  + C Y B 3 ~ Y u ♦ ♦ 8' 

THE  FOLLOWING  ROUTINE  USES  THE  NEWTDN-RfiPHSDN  METHOD  TD  SOLVE  fi 
TRANE  EN DENTAL  EQUATION  SD  THAT  THE  RADIUS  AND  fiNOLE  DF  THE  VERTICAL 
Efib  ENDS  CAN  EE  EfiLCULfiTED. 

THED=3. 14 


554  0 


THE=THED 


64 


IE 


DRflC  IE  ••"01  -'SI 


5550  75  F T H E = B C S V M C ♦ S I N (THE £ . >'  •••" THE- B C S Y M S £ . 

556  0 DFTHE=B  C S Y M C ♦ C □ S (THE  £ . T (£.  ♦THE  7 - B C S Y M C ♦ S IN  (THE  / £ . T T H E ♦ ♦ £ 


cr  c r 

_l . _!  | I I 

553  0 
5590 

56  0 0 
561  0 
5 6 £ 0 
5 6 . ■!  1.1 

564  0 

565  0 

566  0 

567  0 
5680 
569  0 

57  00 
5710 
573  0 

573  0 

574  0 

575  0C 

576  06 
5770 


DELL  TH=-FTHE-"'DF  THE 
THE=THE+DELTH 

I F (THE . GT . £ . *3 . 1 4 1 593 . DP . THE . LT . 0 , 60  TD  516 
IF  CABS  (DELTHT  . LE.  EPSL0N7  SD  TD  76 
60  TD  75 

76  R $ Y M = B C S Y M 6 -r'  THE 

THE  FOLLOWING  STATEMENTS:  ARE  USED  TD  CALCULATE  THE  RIRBRG  VOLUME 
S L □ P E = (.  Y B S - Y C > ••••■  C 7 B S - 7 C 7 

V 1=  SLOP £♦♦£♦  C 7 C ♦ ♦ 3 - ' 3 . - 7 B S ♦ ♦ 3 3 . + 7 B S ♦♦£♦  7 6 - 7 6 ♦ ♦ £ ♦ 7 B S ? 

V £ = Y C ♦♦£♦  ( 7 6 * 7 B S >' 

V 3 =£;  ♦ Y 6 ♦ S LDP  E + ( 7 B S ♦ 7 C - 7 B S ♦ ♦ £ -p'  £ , - 7 C * ♦ £ £ . 7 

V H B 6 D = R B S ( V 1 + V £ + V 3 >'  ♦ 3 . 14159 
ALPHRS=ATAN  C S LOPE 7 

P 6 = Y 6 -*■  P S Y M ♦ 6 D S (THE £ . +RLPHR S 7 

P S = P 6 + 4 . 3 ; ♦ CR S Y M ♦ S I N ( THE £ . 7 ♦ ♦ 3 > ♦ 6 □ S (ALPHA S 7 / (THE-SIN  (THE > 7 

V B 6 = £ . ♦ 3 . 14159*  P S * ( R S Y M * ♦ £ r'  £;  * (THE-1-  S I N (THE  >'  7 7 

V B R I M =3;  14159  ■r  6 . ♦ ( S C + 7 6 7 * ( 3 . ♦ Y 6 ♦♦£+  ( 3 6 + 7 6 V ♦ * £ 7 

V □ L = V fi  E 6 D + V B C + ’• E P I M 

CONFINE  ERG  VOLUME  TD  ORIGINAL  VOLUME  IF  THE  BAG  PRESSURE  IS  LET 
THRU  AMBIENT. 

I F (VOL . 6E . VOLD : AND . PG I . LE . 0 . >'  VOL  = VOLD 


573 


0 


GD  TD  101 


58306  C0MPLE7  ROOTS  (NO  TORSO  AMD  BAG  CONTACT? 

5340  80  VOL = VOLD 

5350  FTT=0. 


586  0 

FTH=0. 

537  06 

60MPUTE  GAS  FLOW  INTO  BAG 

583  0 

1 0 1 GAL  L 6 A S I N ( 7 « 0 1 N 7 

589  06 

SINGE  SUBROUTINE  ° SOLVE*  GALLS  M BRGSUB " TWI6E 

PEP  SDLUTIDN 

59  0 06 

POINT  WE  MUST  DIVIDE  THE  TIME  STEP  BY  £. 

591  0 

DELT  AT  = $TEP-r'£ . 

59  £ 0 

GW  1 =GW+U I N* DEL TAT 

593  06 

TEST  TD  SEE  IF  BAG  DEPLOYED  YET? IS  BAG  PRESS. L 

ATMOS . PRESS 

5940 

IF  (7.  EQ;  0.  DTIME=.  £5 

595  0 

I F (PA . GE . ATMOP . DP. 7. 67 . DT I MET GD  TD  1 07 

596  0 

GW=GW1 

5970 

p Ei = G W ♦ G T Z ♦ U V D L 

5980 

PG 1 =PA-ATMDP 

599  0 

GT =GTZ 

6 0 0 0 

FFT=0. 

6 0 1 0 

FP=  0 . 

6 08  0 

FTT=0. 

6 03  0 

FHEAD=0. 

6 04  0 

FTH=0. 

T' 


65 


- 13  - 


DRHC 


1 O 1 ••-8 1 


605  0 
6 06  0 
6 07  OC 

DTIME-7+STEP 
GO  TO  150 

COMPUTE  NEW  TEMP.  AND  PRESS.  DUE  TO 

NET  GAS 

GAIN  IN 

BAG 

6 08  0 
6 1 i 9 0 
6 1 0 0 
61 1 OC 

1 0 r bT  r’  31  i.  GW  + G r + K- 1 n + gtz  + del  t AT  ' G i i 1 1 
p m IT  M = U ♦ G T 7 ♦ G W 1 
p A 7 * P M U M V 0 L 0 

COMPUTE  MEW  GAS  PRESS . AMD  TEMP ; DUE 

TO  POL 

Y TROPIC 

COMP ; 

6 1 8 O 

6 1 3 0 

614  OC 

P A E'  =«=  ( P M U M V 0 L y ♦ ♦ F P M ■ ■ P A 7 ♦ ♦ -r  F P M - 1 . :■ 
GT8=CGT7+  rPA8.-  PA7  :'  ♦♦  ( (FPN-I ; >VFPN>- 
Ei A G V E M 7 I M G C 0 M P U T AT I ONE  T FI R S T C A L C 

. PRESS 

. RATIO 

ACROS 

6 1 5 0 

6 1 6 OC 

P R 7 = A T M 0 P P A 3 

TEST  FOR  CHOKED  FLOUT  ALSO?  IF  PR7> 1 

? BYPASS 

QE7H . IS 

SET  G 

617  0 

6 1 8 0 

IF  -rpR7.LT.PRS::'  GO  TO  103 
IF  (PR7.  GE . I . > GD  TO  110 

6 1 9 C 

0 L 

1 C 
£( 
3 i; 

4( 


hr 
68 
6 c 
6 c 
6 c 
6 c 
68 
6 c 
6 c 
6 
6 
6 


6 L 

7 C 


6 
6 
6 
6 
6 
6 

64  0 L 
64 1 L 
6 4 c L 

643 

644  i 

645  C 
6 4 e,  r 
647  C 
643  r 
6 4 9 1 

65  0 1 
651 1 
6 5 c i. 

653  ( 

654  r 


FM  1 =VC  1 ♦SORT  {.PR7++  (£  . p N I :<  -PR74.4.  f ,:p^  1 + J pN  I 

GO  TO  1 09 
103  FM I =FM8 

C COMPUTE  EXHAUST  FLOW  AMD  RESIDUAL  GAS  WEIGHT 

I 0 9 0 E 7 H = S 9 RT  ( (778.  ♦P M t -r p M 1 - 1 . >'  ♦ A V ♦PA 8 ♦ F M 1 S Q R T ( U ♦ G T 3 ; 
G 1 1] = i ? i.0 1 — U E 7 H ♦ D t.  L T AT 
GD  TU  111 
110  G ill = GUI 

C COMPUTE  PRESS.  AMD  TEMP.  OF  GAS  AFTER  VENTING 

1 1 1 RATI  L1=GU..-  GW  1 


39  0 

PA=PA3+RAT I □♦♦PNI 

3 0 0 

GT =GT 8+RAT 1 0+*  (PN  1 - 1 

31  OC 

COMPUTE  PRESS.  RATIO 

TO 

DETERMINE  WHETHER  GAS  COMP. 

38  OC 

TIME  THRU  LOOPT  THEM 

SET 

PROPER  POLYTROPIC  E7P0MEMT . 

3 3 0 

p P 6=PA8-'PA5 

34  0 

IF  (PR6. LT. 1 . 0001 > GO 

TO 

118 

35  0 

FPM=PN£ 

36  0 

GO  TO  1 1 3 

3 7 0 

118 

FPN=PN3 

38  OC 

C 0 M P U T E E:  A G P P E S S U R E . 

3 9 0 

1 1 3 

PG1=PA-ATM0P 

IF  THE  BAG  PRESSURE  IS  NEGATIVE?  CALL  IT  0.  FOR  BAG  FORCE  CALCS. 
IF  CPGl . LT. 0 . > PG 1 = 0 . 

IF  THE  TORSO  IS  MOT  IN  CONTACT  WITH  THE  BAG - SKIP  THE  BAG  FORCE 
CALCULATION. 

IF  ■^♦♦£-4.  ♦ A+E>  15  0'  150'  115 

THE  BAG  FORCES  ARE  CALCULATED  IN  THE  ME7T  SEVERAL  STATEMENTS. 

115  E N P H I = P G 1 ♦ S 0 R T ( S A ♦ ♦ £ - W E ♦ ■ ♦ 8 - 4 . > ♦ ( S Q R T ( S C ♦ ♦ 8 ♦ ( 1 - W B ♦ ♦ 8 •••"  ( 4 . ♦ 

csa^+8  y y y y ♦ -r  i - bp  •■•g  i y r? . ♦sort  c s a ♦ ♦ 8 - w b ♦ ♦ 8 4 . +sc+*c:+  -r  i -w  b^+8 
& •••■  r 4 . ♦ s A ♦ ♦ 8 y y y v 

FFT =-£ . ♦ENPH I ♦A.BST ♦COS  ,rBE  TAT 
IF  (BETA.LE. 1.5703>  FFT=0. 

F P = P G 1 ♦ W B ♦ A B S T 

I F ( R B Ab+A Ei S T — P M ♦ I 8 . 1 14  0?  140?  1 3 5 
135  IF (ABST+PBAG. GT. RT0PHM8. > GO  TO  136 
RHEAD=  f ABST+RBAG-RN+ 1 £ . > -•"£ . 


66 


14 


DP  fli: 


IE'.--  U1.-SI 


6 ti  y 
66  0 
661 
6 6-  £ 

663 

664 
6 66 


6 6 3 
6 7 0 
671 


674 


i’ 


b 7y 
68  0 
681 


684 


136 
1 37 


140 

141 


GO  TD  137 

RHEflD=  C CRTOPH-RN ♦ I £ „ >' .-•£ . 

HEADT=WHAC  7 

IF  CEP.  LT.  HEADT-S-' . . flfllt.  MEIfluT . bT . 3.  ♦RHEAD*  WH=£.  ♦AHEAD 
FHt.  AD=E . ♦AHEAD*  '■  * 1 1 H ♦ b 1 + F F T ■-  A B ’S  T .!• 

I F C'y  C6 — Y C7 •'  . LT.  U . > F HE  AD=F HFflD+’L  08  C Y C 6.'1  — Y C7*  ■' 
FTH=TNECK-FHEAD*RHEflD  -"IE. 

FF-FP*  CRN#  1£.  -REflG > ..-'flBS T 
FFT =FFT ♦ (RNM2. -ABAC* -ABST 
GO  TO  141 
FHFflD=0. 

fth=tneck 

FCHEST=FFT+FP 
T R fl  N S T 0 R = - T N E C K 
FTT =-FHEAD^RN^CDS  CY  C6>  -Y  i? 

COMPUTE  THE  BENDING  MOMENT 
MO=FCHEST ♦ T 

COMPUTE  THE  AY  I ALL  FORCE  APPLIED  TD  THE  COLUMN. 

F A C = F C H E S T ♦ C 0 S CTHET  AC- THE  TAT)'  +FHEAD+COS  CTHET  AC  - THE  T AH  :< 
COMPUTE  THE  NORMAL  FORCE  APPLIED  TO  THE  COLUMN. 

F N C = F C H E S T ♦ S I N CTHET A C - T H E T A T V + F H E A D ■ ♦ S I N C THETA C-THETAH > 
F N C OL=F N C + G S 3 E . 1 7 * ! .1 C ♦ S I N CTHETA C > 

COMPUTE  THE  REACTIONS  AT  THE  COLUMN  SUPPORT  POINTS. 

R s C = F N CDL+1 . 5+  CMD+F N C 0 L ♦ L. S C / L F W 


> +TRA  N :S  TOP-RFT  ♦FCHEST 1 E 
APPLIED  TO  THE  COLUMN  AT 


THE  WHEEL  HUB 


R1  B A — i!  R1  N! !_'  ♦ *’  L S'  P — L.  0- 1. 
RBF=  C R :S  C ♦ CLBA-LS:  f 


CLEF-LEA* 

CLEF-LEA.-'' 


686  DC: 

COMPUTE 

687  0 

FRCQL=FC 

688  DC 

COMPUTE 

68  y DC 

RESISTED 

68  0 0 

DELF=FAC 

68 1 DC 

COMPARE 

69  £0 

IF  CSCDLC 

683  D 

i f cveoLC 

684  D 

VCOLCO=  0 

686  D 

146 

VCOLC=VC 

6 8 6 D 

I F CVCOLC 

6 8 7 D 

SCOLC=SC 

688  D 

IF  CSCOLC 

688  D 

GO  TO  14 

7 D D D 

146 

VC  DLC=  D . 

7 D 1 D 

scolc=sc 

7 D £ D 0 

COMPUTE 

7 0 3 D 

147 

7 1=7 I- CS 

704  0 

t 1 =Y 1 - CS 

V + E N C 0 L.  ♦ C L s C - L E F > ;• 

V + F N C 0 L ♦ C L. S C -L  BA >' 

F ACOL =F AC + G S 3 E . 1 7 ♦ W C ♦ C □ S C THETA C >' 

THE  FOLLOWING  ROUTINE  COMPUTES  THE  COLUMN  STROKE. 
STROKE  OF  THE  COLUMN  AND  "VCD LCD"  IS  THE  STROKING 
LOOK  UP  COLUMN  FORCE  AS  A FUNCTION  OF  STROK E . 

CALL  COLF  CSCOL.C  * F COL* 

TOTAL  FORCE  RESISTING  COLUMN  STROKE. 


'S COL CD”  IS 
VELOCITY. 


THE 


DIFFERENCE  BETWEEN  FORCE  APPLIED  TO  COLUMN  AND  THAT  FORCE 
BY  COLUMN. 

IL-F  PC  DI- 
AL TUAL  STROKE  TO  STROKING  LIMIT  OF  COLUMN. 

SE. SLIM* GO  TO  146 
GE.  D. > GO  TD  146 


. LT.  D.  :■  VCOLC=N. 

OLCD+V  CQL  CD^DELTAT  + DELF^  DEL.  TAT 
. L T . S COL  C 0 * S C 0 L.  C = S C 0 L.  C 0 


♦Mr 


NEW  VALUES 


FOR 
O'  * 

o ;y  ♦ 


THE  AIRBAG 
OS C THETA  O' 
IN  CTHET AC* 


COORDINATES 


67 


15 


D FAQ 


1 £••-'  0 1 ••-'8  1 


' n i nr 


JPDATE  OLD  VALUES; 


06  M 

150 

Q 0 L_  Q = V □ L 

0 7 0 

PA5=PA 

MS  0 

L E A = L B A Z - S f . Q L Q 

M 9 M 

I F '"LER  „ L E . 0 „ > RQ 

1 0 0 

L F i 1 1 = L F l.i.l  Z - 5:  f . □ L C 

1 1 M 

veDLeo=^enLe 

1 £ M 

$ COLCG=SCQLQ 

1 3 0 

i.i.lB=ii.if:R6T 

1 4 0 

jilH=!ilHRPT 

I 5 M 

F E L3F= V < 5 > - y ■: 8 > 

1 6 M 

CALL  EPF  I NR  •:  SEN  * 

17  0 

FELK  F=Y  t.  5 > -*-Y  <‘8> 

1 8 0 

C R L L $ P F I N Q ( K P A * 

1 9 M 

F 7 = - <:  s F+F  K N E E + F C 

S 0 0 

N=1 

£ 1 M 

RQ  TO  540 

s 8 0 

5 0 0 

7 ST DP =7 

S 3 M 

N=N+I 

£4  M 

PPINT  5 

£5M 

IF  'rN;  ERn  8>  PFINT 

£ 6 M 

RQ  TQ  540 

£ 7 M 

510 

7STQF-7 

£8  0 

N=N+I 

89  M 

PFINT  5 

3 M M 

IF  (N,  Eip » 8>  PFINT 

3 1 0 

RQ  TQ  540 

3£  M 

518 

7STQP=7 

3 3 0 

N=N+ 1 

3 4 0 

PFINT  5 

3 5 0 

IF  CN; EQ. £> PFINT 

SS  M 

RQ  TQ  540 

37M 

514 

7STQP=7 

3 S 0 

N = N+  1 

3 9 0 

PFINT  5 

4 0 0 

IF  -IN.  EQ.  £>  PFINT h 

41  0 

RQ  TQ  540 

4SM 

516 

7 STOP =7 

4 3 0 

N=N+ 1 

44  0 

PFINT  5 

45  0 

IF  OH.  EQ.  8 > PRINT* 

4 6 0 

RQ  TQ  540 

47  0 

5 1 8 

7STQF-7 

4S0 

N = N+  1 

4 9 0 

PFINT  5 

5 0 0 

IF  <N.  EQ.  87  FEINT  * 

51  0 

& STOPPED . 

58  0 

RQ  TQ  540 

530 

58  0 

7STQF-7 

54  0 

N=N+ 1 

M.  CM  SF  « NSF 


1 7 V 


68 


16 


dp  He  lp.-i-ji  r8i 

755  ft  PR  I NT  5 

756  0 IF  < N . E 0 .•  8 Y P P I N T * " p H I N □ T p □ N V E R 6 1 N 6 1 P U N 5:  T 0 P PE  ft:  “ 

757 ft  54  ft  RETURN 

758  Fi  END 

759  OP 

F1  -•  D 0 6 THI  5 8 IJ B 8' □ u T I N E C 0 f’l P i.:'l T F $ 7 HE  PR T F T H R T G R S'  F N T F p x THE  J: R G • 

761  ft  S U B R □ U T I N F G h 8 I N f 7 « ijj  I n 

7 6 £ 0 e 0 M M o N ■r 6 ft S FLO -r' N P 6 * ft g N ( 8 , 8 4 

!*'  68  ft  L UL  L L [DDF  UP  1 i HE  N ’ NE’ H i h.i  J N > 

7640  RETURN 

765  ft  E N D 

766  ftC 

767 ft 0 THI S S U B P □ U T I N E 6 □ M P i:l T E S T H E p 0 p U M N F f] R p F ft 5 ft  F p N Q F S T R G K E . 

7 6 8 ft  S G B R G U 7 1 N F G G L F <"  5:  p 0 L C t F P D L ;■ 

7 6 9 ft  e □ M M G N e Q L.  F G R 8 N P p < p 0 L ’■  £ t £ 4 :> 

7 7 ft  ft  p ft  L L L □ G K U P ( S C □ L p « p G L * N P p i F P D L 7 

7710  RETURN 

7780  END 

778  OP 

774  ftp  . THIS  SUBROUTINE!  PL.  ft  PE'S  PERTft  I N VRLUE'S  IN  MRTR  1 7 FORNftT  FOR  PR  I NT  INC 

775  ft  S' U B P 0 U T I N E P PINT!  f 7 < Y ■»  0 Y > 

776  ft  p: □ M M 0 N 8 [] U T 1 8 N6  * 7 6 '■  4 1 5 ft 7 « T 6 ( 5 ft  :> 

777  ft  p 0 M M □ N 8 H I p 8 T H I P 1 7 5 >'  * H P 6 S > 1 7 5 > < p I N 7 1 

7 7 8 ft  D 0 U B L E P R F P I P 1 0 N f ( £ Y 

779 ft  D I M F N S I □ N D t f 8 > 

7 8 ft  ft  p G M M 0 N FI H N D ft T 8 7 p « 7 T <•  7 H * R T * R N * RH  h RT G P H h X £ Z * Y 8 2 ii  Ml B 

781ft  N 6 = N E,  + 1 

7 8 8 ft  I F ■:  N 6 „ ft  T „ ft  ft  Y p E T U E:  N 

7 8 3 ft  P H = P G S' ,r  Y ( 6 '•  1 

784ft  PT =POS  ,rY  <■?  :.< 

7 8 5 ft  S H = S I N ,r  Y C 6 7 7 

7 8 6 ft  S T = S I N ('  Y (.  7 7'  1 

787ft  76  f 1 n N6  ■■  = CDY  ( 1 > ♦PT +PT+DY  (£.7  :<  88.9  D l 7 

7 8 8 0 7 6 C 8 * N 6 7 = ( D Y ( 1 7 ♦ S T-*-RT ♦ Y ( 3 7 ♦ Y C 3 7 7 •r  3 8:  17 

7 8 9 ft  7 6 ( 3 N 6 1 = RH+D Y F 8 7 + D Y 1 1 7 * C H+R N ♦ ( D Y 3 •'  * '■  p El ♦ P: T + x. E! x 1 V + Y C 3 > 

7 9 ft  ft  ;?■  ♦ t 'r  3 7 * <’  S H ♦ p T - p H ♦ S:  T 7 7 7 8 8 . 17 

791  ft  76  ,r4 1 N6>'  = (PH^-Y  C.£7  ♦Y  fS.-1’  +RN*  Pf  '13  * f (37  * PSH^ST+PH^GTT  — D'Y  P8  -1 

7 9 8 ft  ♦ < S El  ♦ p T - p H ♦ S T Y > - [i  Y 1 > ♦ S H Y 8 3 8 . 17 

7 9 3 ft  T 6 t N 6 > = 7 ♦ 1 ft  ft  ft : 

794 ft  THI P f N6 > = (N6*P I NT  1 Y ♦♦ ft . 6 

7 9 5 0 H P G S v N 6 '■  = S Q R T ( 7 6 C 3 * N 6 > ♦♦8+ 7 6 4 1 N 6 ♦♦  8 7 

796ft  RETURN 

797ft  END 

798  ftp 

799  ftp  THIS  SUBROUTINE  SOLVES  THE  DIFFERENT  I RL  FOUNT I ONS  THRT  DETERMINE 

8ft  ft  ftp  THE  DRIVER  KINEMPTIPS  RND  VEHIPLE  MOTION:  THE  FOURTH  ORDER  RLlNGE- 

8 ft  1 ft  C , K IJTTR  M FT  HOD  I S U S E D T 0 S TRRT  THE!  I N TEGRRT 1 0 N " B U T 0 N p E THE  FIRST 

8 ft 8 OP  FOUR  POINTS  RRE  OBTAINED  M1E  SUITPH  TO  THE  MORE  EPDNOMIPHL  FOURTH 

8 03 ftp  ORDER  RDRMS ^MOULTON  PREP I PTOR-PORREPTOR  METHOD. 

8 ft  4 ft  S U B P 0 U T I N E S 0 L.  V E ( N 


69 


- 17 


HP  HU 

1 8‘7  0 1 .--8 1 

8 05  0 

C □ M N □ N 7 TINE 7 S TEP  < 7 S TOP 

806  0 

UDNNGN 7(4 1 U 7 T H I C U I 75?  * HPRS  C 1 7 

8 07  0 

D □ U BLE  PPE U I S I □ N Y •: 8 ? n 't  T 1 8 ? 

ft  08  0 

PERU  Ot 

8 0 8 0 

0 □ U B L E P P E U I S I □ N B £ 7 0 1 B 1 8 n B £ 

SI  0 0 

0 1 N E N S I □ N 0 t -r  8 ? - P -r  3 , 8 ? , . 1 f 8 ? , 

3 1 i 0 

B £ 5 1 =85 1 . 

8 1 ,7  0 

6870=870. 

S I 3 0 

618=18.- 

5: 1 4 0 

PPI 1=0. 

3 1 5 0 

PP I £=  0 . 

8 1 P 0 

73T0P=. 140 

8 1 7 0 

STEP= . 001 

8 1 8 0 

URLL  SETUP  CYi  Y ? 

8 1 '9  0 

0 R L L 0 1 E E U C 7 * ?'  * 0 Y :■ 

88  0 0 

P I NT  1 = . 0 05 

88 1 0 

P I NT  8= . 0 05 

888  0 

URLL  PP I NT  1 Of , Y i OY? 

883  0 

URLL  PRINTS  C7-  Y- DY> 

884  0 

I F CYST OP . RT . . 85? 7ST0P= . 85 

385  OU 

STRPT  OF  I NTE8PRT I ON  ROUTINE 

88  P Oh 

PUNSE  KUTTR  STRPT  UP 

887  0 

J 1 ? = 1 

88  8 0 

J US? =8 

383  0 

J C3>  =3 

8 3 0 0 

DO  787  0 kr=lnN 

8 3 1 0 

~P  I'  f | 
1 L_  1 U 

F ( 3 h 7 ? =DY  C7? 

838  0 

00  7550  J7= 1 h 3 

833  0 

00  7 3 0 0 7=1 < N 

834  0 

7 3 0 0 

S C7?=0Y  CO  ♦STEP 

835  0 

:yn=:y+step7£; 

S3P0 

DO  733  0 7=1  n N 

837  0 

733  0 

YT  ,r7 ? =Y  (K?  +S  < 7?  78; 

838  0 

U R L L D I F E O 'r 7 N * YT  - 0 Y ? 

838  0 

00  7 3P0  7=1 « N 

84  0 0 

73P  0 

S '.rK>  =S  ,r7?  +8.  ♦OY  CK?  ♦STE  P 

841  0 

00  733  0 7=1 h N 

848  0 

738  0 

Y T <■  7 ? = Y '■  7 ? + S T E P ♦ 0 Y < 7 ? 7 8 . 

843  0 

U RLL  01  FEU!  CYN < YT n OY? 

844  0 

00  7410  7=1 n N 

845  0 

741  0 

3 y 7 ? = S ,r 7 ? + 8 . ♦ 0 Y <■  7 ? ♦STE P 

8 4 P 0 

00  7430  7=1 « N 

8470 

74  30 

YT  (7?  =Y  -17?  + DY  U7?  ♦STEP 

848  0 

7=7  + STEEv 

849  0 

URLL  01  FEU  CYh  YT n OY? 

85  0 0 

00  747  0 7=1 « N 

851  0 

747  0 

■f  r f v = y •:  7 ? + f s y 7 ? + 0 Y y 7 ? ♦ S T E p ? -r 

858  0 

URLL  01  FEU  U7? Y- OY? 

85  3 0 

ROTO C 7 5 0 On  75 3 0-75 5 0 > - J 7 

854  0 

75  0 0 

00  7510  7=1 < N 

' * PI  NT  I 
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r 


- is  - 


hr  fie 


1 £-'•  O 1 .-s  i 


8550 

7510 

F (2*  K?=DY  (K>‘ 

856  0 

60  TO  7550 

8570 

7530 

HO  7540  K=1#N 

85  SO 

754  0 

F(OK)=Ht(K)' 

85  SO 

755  0 

RRI1=Y 

86  0 0 

PR I 2=7 

861  OC 

P R E B I C TOR-CO R R E C T 0 R S E C T 1 0 N 

868  Oe 

F'REH  I CTOR 

863  0 

759  0 

HD  76  00  K=1 « H 

S64  0 

760  0 

t T (K>  = t CK>  +STEP*<55;  ♦BY  OO'  ^5 

865  0 

&-9.  ♦F  (.J  *r 3 > t K >'  > -r'2'4 . 

S66  OC 

SAVE  H’f  •*'  S 

S67  0 

HO  764  0 K= 1 ? H 

86  S 0 

764  0 

F ( J C 3 :•  , K>‘  = D f ( K>' 

869  OC 

EVALUATE  STEP 

870  0 

f=Y+STEP 

S71  0 

CALL  HI FEC!  (ft  YT*  HYY 

87 £ OC 

ROTATE  VECTOR  RD INTER 

873  0 

IJT = J (3>‘ 

S74  0 

j ox = j .r £ ;y 

875  0 

j C2>'  =j  i >• 

S76  0 

j •:  i y =ut 

S77  OC 

CORRECTOR 

S7S  0 

HO  7750  K= 1 t H 

8790 

775  0 

f ( K >•  = f ( K >■  + S T E P ♦ ( 9 . ♦ H Y ( K >'  + 19. 

88  0 0 

&K>  y .-'24 . 

SS£  OC 

ADDITION  OF  ERROR  TERM 

S3  3 0 

HO  78  0 0 K=1 « N 

SS4  0 

780  0 

f ao'  = CB25 1 ♦f  aa +E 1 9+YT  (K?  y .-B 

885  OC 

SECOMH  EVALUATION  STEP 

836  0 

CALL  H I FEU! ( f t t t H Y •' 

SS65 

CALL  UP  HATE  (ft  Y t DY7 

SS7  OC 

PR I NT I NS  SECTION 

S3  8 0 

PR I 1 =PR I 1 +STER 

3890 

PR I 2=PR I 2+STER 

39  0 0 

I F (PR 1 1 . LT . R I NT  1 >' 60T07S9 0 

39 1 0 

PR  1 1 =PR  1 1 -F'  I NT  1 

89  £0 

CALL  PRINT  1 (fit  YnHYV 

393  0 

789  0 

I F (PR  12;  LT:  RI NT2>' 60T0792 0 

894  0 

PR  1 2=PR  1 2--P  I NT  2 

3950 

CALL  PR  I NT 2 ( Y * Y t BY? 

896  0 

7920 

IF  (f. LT. YSTOP> 60T07590 

3970 

8 0 0 0 

CALL  PR  INTI  (fit  ft  BY)' 

393  0 

CALL  PRINT2 (ft  ft  HYT 

399  fi 

RETURN 

9 O 0 0 

ENH 

9 0 1 OC 

9 02  OC 

THIS  SUBROUTINE  COMPUTES  THE 

9030C 

SEAT  FRICTION  FORCE.  HYSTERE 

904  0 

SUEROUT I NE  SPR I NG (F , BELT A - H I 

TRAINT  CRUSH  FORCE  fiMD 
■ TS  CfiN  EE  INCLUDED. 

' t SLOPES  * FORCE*  NPTS> 


THE 


71 


19 


DRRC 

iTi 

y 

OJ 

T— ! 

905  0' 

DIMENSION  F (Eu 84  > 

90  SO 

R=0. 

9 0 7 0 

I F C D I $ T . 6 £ . D E l:.  T R > 6 0 T □ 8 3 

9 08  0 

M = S 

9 0 9 0 

UT=D 1ST 

9 1 0 0 

G □ TO  8360 

9 1 1 0 

8 1 1 0 

M = 3 

9 1 8 o 

UT  =BELTR 

9 1 3 0 

6 [ j T □ 8 3 6. 1. 1 

9 1 4 0 

8 1 4 0 

F 1 1 K I = DEL  TP 

915  0 

F C 8 ;•  K > = F □ R C E 

9 1 S 0 

IF  fk  . i5 T.  5' 6 □ T □ 8 1 8 0 

9 1 7 0 

i5  0 T 0 f 8 8 4 0 * 8 8 4 0 < 8 3 0 0 j « k 

9 1 8 0 

8 1 8 0 

k k =k -3 

9 1 9 0 

R=  1 „ 

98  0 0 

DO  8810  L=8nNFTS 

98 1 0 

F 1 i L ■'  =F  1 < L+k  k 

988  0 

881  0 

F f 8 * L ■'  = F '■  8 * L:  + k ) 

988  0 

N FT  S=NRTS--k  k 

9 .9  4 0 

GO  TO  8300 

985  0 

884  0 

k k =3-k 

9 8 6 0 

R=8. 

987  0 

DO  8880  LL=1hNFTS 

988  0 

L=NFTS+ 1 -LL 

989  0 

Fvl-i  L+KIO  =6  irN  L :• 

9 8 0 0 

8 8 8 0 

F <■'  8 « L + K K ■'  = F f£«  L •' 

9 3 1 0 

NRT3=NPTS+k'K 

938  0 

8 9 o 0 

F f i * 8 ) = F II  * 5 ■'  — F 1 8 * 3 •'  S li  0 R 

9 8 3 0 

F '39 1 8 :•  =0. 

934  0 

F (£<  I >'  =-SLQPE8+F  H 1 8 > 

935  0 

f f 1 h i y = o . 

936  0 

854  0 

U=1 

9 3 7 0 

1J  T = D I S T 

938  0 

836  0 

DO  8370  4=2 * NFT3 

9 3 9 0 

837  0 

IF  IF  ft « jy  . GT.  LIT >' GO TO 8 380 

94  0 0 

8 58  0 

K=  :.F  I 

941  0 

8 39  0 

F 0 R C E=  f 0 T j-F'  f 1 * K y y ♦ f F f 8 ■*  J 1 

948  0 

IF  fR. EO.  1 . >' NRTS=NPTS+KK 

94  3 0 

IFfR.E © . ,3 . N P T 5: = N P T 9: - K K 

944  0 

G □ T 0 f 8 4 0 -•  84 1 0 t8I4 0 < M 

945  0 

841  0 

IF  IF  □ R C E . L E . 0 . 1 G 0 T □ 8 4 5 0 

946  0 

IF  f w B S f f F f£«  J y — F f 8 * K >'  V f F 

947  0 

CG0T0845 n 

948  0 

GO  TO  8110 

9 4 9 0 

.'3  4 5 0 

RETURN 

95  00 

END 

95 1 OC 

958  OC 

THIS  SUBROUTINE  FORCES  CE 

95:9  0 

S 1 J B R 0 U T I N E F R I N T £ f 7 « Y < D Y ? 

954  0 

REAL  MO 

PR  I NT  I N6 


72 


f 

f 

) 


£0  - 


HRRl: 

1 £■■"  01. ••-81 

9 5 5 0 

COMMON  •••'□UTVN 4 T 735  > - '?■  0 C6»  35  T , 71  76.4  35 > 1 7£  C6? 

3 5 7 4 7 3 7 6 < 3 5 7 4 7 4 7 F 

9 5 6 0 

75  '■  5 1 357 

9 5 7 0 

DOUBLE  PRECISION  Y 78  • 

9 5 9 0 

DIMENSION  D’f  797 

9 5 9 0 

COMMON  ■•'T’lwN DPT  -r'7L  4 £T  4 £N  <1  PT  4 PN  4 PH  4 PT  OPH  4 ’7 £7  4 

Y£Z 4 MB 

9 6 0 0 

COMMON  r'pHPRM 'PpT  * THE  ThT  Si  THEThHZ 

96 1 0 

COMMON  *MPhPhM..-FhCOLi  S COL Ci  PG1 « FKNEE9  SFY  PNC 

OL 4 MO" PPCOL * VCQLC 

96  £ 0 

C 0 M M 0 N M M pRPfi M E P , V 0 L - F FT  - F p 

96  3 0 

D H T w F 1 3 • r h r „ 6 8 1 8 1 8 1 8 • - G ••-' 3 £.17 * D ••■- 57.295 

7£0.-- 

964  0 

N = N + I 

9 6 7 0 

IF  fN:  FT.  35':'  PE  TOP  N 

96  f f i 

CT=C'0  7 Cf  7 7 :•  :• 

9 7 0 

T 7 N T = 7 ♦ 1 0 0 0 „ 

9 6 £ 0 

17  0 7 I 4 N t 74  :>  .--'G 

9 6 Ff 

:7  0 '79  h N:<  =Y  C 4 ■ 

9 7 0 u 

7 0 03'  NT  = V 7 9 :>  +F 

971  0 

% m i 4- 1 N * = ( 1 ! y ,r  1 j y T + PT ♦ [!  V '*  1 ' y 

97  £U 

7 0 '"5'  N :*  = SC DEC 

97  3 0 

7 0 *16 1 N t =P0 1 

9 74  0 

7 1 7 1 1 N ■'  = '■(  7 5 ♦ F 

A t ' "j  M 

7 1 7 £ ' N 1 = Y 7 1 7 ♦ v 

9 7 6 0 

7 1 7 7;  h N ■>  = — D 7 7 1 7 ,r'  0 

977  0 

7 1 74'  N :•  = E K M F E ■r'  £ . 

978  0 

7 I 7 5 ' N 1 = S F 

9 7 9 0 

7 1 7 6,  4 N 1 = if  V 7 5 '>  7 >r  8 T 7 ♦ F 

9 ? 0 0 

7£  7 I 4 N 7 — Cf  (5  ! +PT ♦ OS  I N 7 f 777  7 -*-S  I N 7 THE  T 97  ,77  :>  :< 

♦ F 

9 8 I 0 

7 £ 7 £ ' N ■>  = f 7 7 7 ♦ D 

9 3 £ 0 

7 £ 7 3 « N 1 = Y 7 3 7 ♦ D 

Q '•**!  y,  r 1 

7 £ 7 4 ' N :>  = Ei  Y 7 3 7 ♦ D 

a y *4-  m 

7 £ 7 5 4 N 7 = 7 £ 7 1 4 N •+■  Y 7 3 1 ♦ F 

9 6 5 0 

’•m  y h-.  f-"j 

7£  76  < N :<  = 7PT ♦ Y 7 :3  , ♦FT  + 7’  7 1 < j-‘f  74  .1  7 ♦*.;? 

,7  3 7 1 4 N * = F ♦ 7 Y 7 5 7 + P N ♦ 7 S I N 7 Y 7 7 7 7’  S'  I N 7 T H F T 8 T 7 7 

7 +PH+  7:3  J N 7Y  76  :■  7 — 

9 q 7 1 j 

7 3 I N 7 THE  T 9 HZ 7 7'  7 

'9  7:  8 0 

,7  3 7 £ 4 N 1 = V’  7 E,  :•  ♦ 

9 9 q |*| 

‘73  78'  N 7 =Y  7 £7  ♦[! 

9 9 I j | j 

;7  3 7 4 ' N 7 = D Y 7 3 ;i  ♦ f_i 

9 9 1 0 

■7:3  75  4 N 1 = 7:3  7 t 4 NT  -*-'f  7 3 > ♦£ 

y :a  y m 

:7 3 7 E,  4 N T = 7 Y 7 6 7 — Y 77 7 7 ♦ D 

y y f-'i 

74  71  4 NT  =FhCDL. 

9 9 4 0 

74  7£ 4 N ■ =FNCOL 

9 97  0 

74  7 3 4 NT  =M0 

A A y || 

74 74' NT =FPCOL 

997  0 

74  754  NT  = SC DL.C 

9 9 8 1 J 

y 4*  f t*  < f‘T  1 = i,m?  i.i’  n Ll  L 

9 9 9 0 

75  7 1 4 N . =BP 

1 0 0 0 0 

75 72 < N 7 = VOL 

1 0 0 1 0 

75  73 « N >'  =901 

1 0 0,9  0 

75  7 4 h NT  =FFT 

1 0 0 3 0 

75  75  4 N T =FP 

1 0 0 4 0 

PETOPN 

1 0 05  0 

END 

1 0 06  0C 
1 0 07  0 

S 1 J B R 0 U T I N E U F'DRTE  < 7 4 Y 4 D Y 7 
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DR  HU 


1 0 080 

DOUBLE  PRECISION  Y <8> 

1 0 09  0 

DIMENSION  DV <8> 

1 i J 1 0 0 

C D M M 0 N - K N E E R E S T - N K R < S K P * R U T ' REL K R 

1011  0 

C 0 M M 0 N SERT F P'  I C •-  N S F « S U N ■.  S PUT  ? RE LSF 

1 0 1 £ 0 

SFUT =REL.S'F 

1 0 1 3 0 

RUT =RELKR 

1014  0 

RETURN 

1015  0 

END 
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Appendix  E - DRAC  Sample  Run 


Ei  R A C 0 8:41  fc.  $ T 1 8 0 £ 8 1 

INPUT  FILE  NAME^LRSV 


INPUT  URL UEJ 


I N P id  T U N I T S '■  N 3 E C $ r-1  P H * fi E 6 R E E $ i I N L A E 3' * LBS  ^ FT -‘•LB -3  * 6 J 3 


INITIAL  VELOCITY: 

INI TIRL  HERD  R NEELE" : 
INITIAL  TORSO  ANGLE: 


:>  „ 0 0 
-18.  0 


ML  EG 

M TOR 30  MHEAD 

F'T 

RN 

RH 

RTDPH 

77.  0 

67.0  11. 0 

1 3 „ 8 

18.5 

4.  0 0 

88.  0 

NETS  NEC 

K NPTS  KR  NETS  UEH 

NETS  S 

E A T N P 7 S G A S 

NPTS  C 

□L  SL.ST 

SL . KR 

? 

0 fit. 

Ei 

£ 

0 . 1 0 0E+  04  0 . 

88  OE + 04 

1 ^ fl  S h 

LOU  TIME 

- 1 0 0 . 

14. 0 16. 0 

1 8 . 0 

38.  0 

bO  ■ U 

88.  0 

1 0 0 . 

GAS  F 

LOW  - LE 

y v p C' 

0 . 

ij.  4. 88 

' -!  o j-jK, 

3.  07 

0 . E’  I 4 

0 . 

0. 

C OLL'IN 

N STROKE 

- INCHES 

- 1 0 0 » 

0 . 0 . 8 5 1 j 

Ij  . Li  Ij  Ij 

I . 5 1 j 

8.  0 0 

COLUMN  FORCE 

- LBS 

0 . 

0 . 4 0 1 j „ 

4 0 0 . 

0.160E+04  0. 

1 E OE  + 04 

SEAT 

FRICTION 

h j i p [_  H 1.;  p U p H i 

- 1 0 0 . 

0 . 0 , 5 0 0 

I 5 . 0 

1 6 . 0 

1 0 0 , 

SEAT 

FRICTION 

FORCE  - L.B3 

0 „ 

0 . 4 o 0 . 

3 0 0 „ 

0 . 

0 . 

NECK 

ANGL  E 

-8  0.  0 

-75. 0 -60. 0 

-3  0.  0 

0 . 

3 0.  0 

60.  0 

75.  0 

8 o . 0 

NECK 

TORCUE 

8 0 0 „ 

15  0.  1 0 0 . 

5 0 . 0 

0 . 

-50.  0 

- 1 0 0 . 

- 1 Li  0 . 

-8  0 0 . 

UEH. 

p U L 3-  E — 

T I ME 

- 1 0 0 . 

0.  40.  0 

I 08. 

115. 

1 5 0 . 

UEH. 

P ij  o.  S E — 

DECELERAT I QN 

0 . 

0 . 81. 0 

■?i  :<  . cr 

L_  .Jt.  D _l 

0 . 

0 . 

KNEE 

D I 3 P L A C E 

ME  NT 

-5  0.  0 

8.00  8.75 

6 „ Li  0 

7 . 5 0 

8.  0 0 

11.0 

50.  0 

KNEE 

FORCE  - 

LBS 

0. 

0 . 4 ij  0 . ij . 

c y Mfc  + IJ4 

0.8806+04  0. 

1 8 OE  + 04 

0 . 1 9 0E+  04  ij . 

1 0 0 E + 0 4 

AT  MOP 

AGS  AT Z 

u 

PNI 

RMS 

PN3 

14.7 

-14.7  0.116E+04 

E,  E'  0 . 

1.40 

1.40 

1.40 

VC  I 

UC8  AU 

SA 

SC 

717 

T 1 Z 

DON 

ij . 7 ij  i j 

0 . 7 0 ij  3 . 5 ij 

11.5 

7 . 5 0 

88.  5 

83.  0 

8.50 

SLIM 

T HE T AC  MU 

}■  vjT-T 

L F ill  Z 

LBAZ 

LBFZ 

ill  IS 

8.  0 0 

17. 0 0 . 1 6 0 

1 6 . 0 

14.  7 

1 6 . 0 

lb  b a C 

1 8 . 0 

l.i.l  H 

'TO  R I MR  AD 

cl  b. 

T q'  .it 

MB 

8.  0 0 

U ■ r m »’  5 

36. 5 

8.  0 0 

17.  0 

75 


INITIft L V ft L ij E S: - C □ N V E R T E D U N I T S < S E C 9 F T S E C - R ft D I ft N S » FT n LBS*  FT- L B S * F T S E C ♦ ♦ £ > 

INITIAL  V Eli  DC  I T'f : 57.6 

I N ITIftL  HE ftB  ft N 6 LEs  - 0 . 8 7 3 E - 0 1 

I HI  T I ft  L T □ R S □ ft  N b L E ■“  ~ 0 . 3 1 4 


ML  ft  ft 

M TOR SO  MHEftD 

RT 

RN 

PH 

RTOPH 

i*1  n , “J  Q 

3 . IJ  ft  IJ » 3 4 c 

1.15 

1 . 63 

0.  333 

3.  33 

MPT $ NEC 

K NETS  KR  NETS  VEH 

MPT"  S 

EftT  NPTS  6 

ft S NPTS  C 

OL  SL.ST 

SL . KR 

9 

f-, 

K 

ft 

6 

0.  130E+05  0. 

34SE+05 

ftft 3 FLQiiil  TIME 

- 0 . I 0 0 0 ■ 

I 4 0 E - 0 1 0 . 1 6 0 E - 0 1 0 . 

1 ft  OE  - 0 1 

0. 33  OE- 0 1 

0 . 65  OE- 0 1 

0 . ft 3 OE  - 0 1 

0 . I 0 0 

ft  ft  3 FLOM  - LB 

•r'  S E C 

0. 

0.  4. ftft 

3 a 86 

3.  07 

0 .614 

0. 

0. 

COLUMN  STROKE 

- INCHES: 

- 1 ij  ij  „ 

0.  0.350 

0 . 5 0 0 

1 . 5 0 

ft.  00 

COLUMN  FORCE 

- LBS 

0. 

0 . 4 fj  fj  „ 

4 0 0 . 

0. 160E+04 

0 . 160E+04 

* 

SEftT  FRICTION 

D I S'  P !.:  ft  C E M E M T 

-ft  „•  33 

0.  0.41 7E-0I 

1 . 35 

1 . 33 

|ij  ■ «.J ' tj 

SEftT  FRICTION 

FORCE  - LBS 

0.- 

0 . 4 0 0 „ 

3 0 0 . 

0. 

0. 

NECK  ftNftLE 

-1 . 57 

-1.3 1 -1.05 

-0. 534 

0 . 

0.  534 

1 . 05 

1.31 

1 . 57 

NECK  TOR COE 

£'  Hi  0 . 

1 5 0 . I 0 0 . 

5 0.  0 

0 . 

-50.  0 

- 1 0 0 . 

- 1 5 0 . 

-*-3  0 0 n 

v'EH.  PULSE  - 

T I ME 

- ij  „ I ij  U 
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